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Abstract 


Fine-  and  microstructure  profiles  collected  over  Fieberling  Seamount 
at  32®26'N  in  the  eastern  North  Pacific  reveal  three  intensified  phenomena 
co-existing  in  a  200-m  thick  layer  atop  the  summit  plain:  (i)  an  anticyclonic 
vortex  cap  of  core  relative  vorticity  -0.50/,  (ii)  diurnal  shear  of 
15  cm  S“V40  m,  and  (iii)  turbulence  levels  corresponding  to  an  eddy 
diffusivity  ATp  ^  10  X  10-4  n;i2  s-1.  The  vortex  cannot  be  explained  solely  by 
Taylor-Proudman  dynamics  because  of  its  -0.3/V2  negative  potential 
vorticity  anomaly.  The  fortnightly  cycle  in  the  vortex's  strength  suggests 
that  it  is  at  least  partially  maintained  against  dissipative  erosion  by 
rectification  of  diurnal  tides.  Baroclinic  diurnal  fluctuations  are  slightly 
subinertial  and  the  diurnal  kinetic  energy  also  tracks  the  beating  of  the 
barotropic  Ki  and  Oi  diurnal  tides.  Their  horizontal  stmcture  closely 
resembles  a  seamount-trapped  topographic  wave.  However, 
counterclockwise  turning  with  depth  of  the  horizontal  velocity  vector  and  the 
180°  phase  difference  between  radial  velocity  u/  and  vertical  displacement 

=  -T/Tz  (producing  a  net  positive  radial  heat-flux  <UrT'>)  are  more 

consistent  with  an  upward-propagating  vortex-trapped  near-inertial  internal 
wave.  The  slightly  subinertial  frequency  is  allowed  by  the  strong  negative 
vorticity  of  the  vortex  cap.  A  vortex-trapped  wave  would  encounter  a 
vertical  critical  layer  at  the  top  of  the  cap  where  its  energy  would  be  lost  to 
turbulence.  High  observed  turbulence  levels  of  e  =  3  X  10-8  w  kg-1  imply 
decay  times  for  the  wave  and  vortex  of  less  than  3  days,  emphasizing  the 
strongly  forced  nature  of  the  system.  Inferred  eddy  diffusivities  two  orders 
of  magnitude  Lrger  than  those  found  in  the  ocean  interior  suggest  that, 
locally,  seamounts  are  important  sites  for  diapycnal  transport.  On  basin 
scales,  however,  there  are  too  few  seamounts  penetrating  the  main 
pycnocline  to  produce  a  basin-averaged  diffusivity  of  0(10^  m2  s-1). 

Interior  and  near-bottom  mixing  likely  contribute  comparable  diapycnal 
buoyancy-fluxes  at  these  depths. 


1.  lntroduction 


The  peaks  of  seamounts  have  been  found  to  be  sites  of  enhanced  near¬ 
bottom  subinertial,  diurnal,  inertial  and  semidiurnal  currents  (Noble  et  al 
1988;  Noble  and  MuUineaux  1989;  Genin  et  al.  1989;  Padman  et  al.  1992; 
Eriksen  1991;  Kunze  and  Sanford  1986;  Kunze  et  al.  1992).  These  elevated 
velocities  support  abundant  populations  of  benthic  filterfeeders  (Genin  et  al. 
1986;  1992;  Levin  et  al.  1994).  How  topography  intensifies  these  currents  is 
is  not  established.  Kunze  and  Sanford  (1986)  interpreted  their  signal  as  a 
critically-reflected  upward-propagating  near-inertial  wave  based  on  the 
counterclockwise  turning  of  its  horizontal  velocity  with  depth.  On  the  other 
hand,  Genin  et  al.  (1989),  Eriksen  (1991 ;  1995b),  Noble  et  al.  (1994)  and 
Brink  (1995a)  report  slightly  subinertial  diurnal  tidal  frequencies  associated 
with  the  dominant  motions  above  Fieberling  Guyot  so  they  advanced  that  the 
fluctuations  might  be  seamount-trapped  topographic  waves  (Brink  1989; 
1990). 

Sloping  topography  allows  bottom-trapped  topographic  waves  with 
frequencies  co  <  A^sina  (Rhines  1970;  Huthnance  1978).  Sloping  bottom 
topography  will  also  cause  internal  waves  with  characteristic  slopes 
CgzICgH  =  knfkz  identical  to  the  bottom  slope  a  to  be  critically  reflected 
(Wunsch  1969;  Phillips  1977;  Eriksen  1982;  1985)  and  scattered  to  high 
wavenumbers  (Baines  1971;  BeU  1975;  Gilbert  and  Garrett  1989;  Muller  and 
Xu  1992),  dramatically  amplifying  internal- wave  shear  and  strain. 

Enhanced  fmescale  shear  and  strain  wiU  support  stronger  and  more  frequent 
shear/advective  instability  to  produce  intensified  turbulence  (Henyey  et  al. 
1986;  Gregg  1989;  PoLzin  et  al.  1995).  Relative  to  the  ocean  interior, 
enhanced  turbulence  has  been  documented  over  seamounts  (Nabatov  and 
Ozmidov  1988;  Lueck  and  Mudge  1996;  Toole  et  al.  1996).  This  has  led  to 
speculation  that  seamounts  act  as  stirring  rods  for  the  world  ocean,  with 
topographically-induced  mixing  near  boundaries  dominating  over  the  weak 
diapycnal  mixing  of  the  ocean  interior  (Lueck  and  Mudge  1996). 

Mean  circulations  are  also  expected  over  seamounts.  Both  Taylor- 
Proudman  dynamics  (Hogg  1973;  Swaters  and  Mysak  1985;  Roden  1987)  and 
topographic  rectification  processes  (Loder  1980;  Maas  and  Zimmerman 


1989a;  1989b;  Haidvogel  et  al  1993)  could  lead  to  formation  of  a  vortex  cap 
over  a  seamount's  summit.  Stratification  limits  the  vertical  extent  of  vortices 
above  seamounts  (Zyryanov  1981;  Chapman  and  Haidvogel  1992)  to 
H  =fLIN.  For  Coriolis  frequency/=  7.8  X  10-5  s-1,  buoyancy  frequency 
N  =  43x  10-3  s-1  and  seamount  radius  L-1  km,  the  cap  thickness 
//~100m.  A  Taylor  cap  would  be  anticyclonic  and  requires  an  impinging 
geostrophic  flow.  In  the  absence  of  damping,  this  vortex  would  have  the 
same  potential  vorticity  as  the  surrounding  ocean.  The  properties  of  a 
rectification-driven  vortex  depend  on  the  nature  of  the  forcing  and  damping 
and  are  less  certain.  They  need  neither  be  anticyclonic  nor  preserve  potential 
vorticity. 

To  better  understand  the  impact  of  seamounts  on  internal  waves,  tides 
and  turbulent  mixing,  fine-  and  microstructure  profiles  were  collected  over 
Fieberling  Guyot  during  March  1991  (Montgomery  and  Toole  1993).  Here 
we  use  the  profile  time-series  and  surveys  to  characterize  the  temporal  and 
spatial  structure  of  motions  atop  Fieberling  Guyot’s  summit  plain.  Toole  et 
al.  (1996)  examine  the  turbulent  boundary  layer  above  the  seamount’s  flanks. 
We  fmd  a  vortex  cap  atop  Fieberling  Guyot  that  is  ~  2(X)-m  thick,  has  a  core 
vorticity  C  =  -  (0.50  ±  0.1 5)/ and  a  core  potential  vorticity  anomaly  of 

The  presence  of  this  anticyclonic  vortex  raises  the  possibility  that 
the  diurnal  oscillations  coinciding  with  the  vortex  are  not  seamount-trapped 
topographic  waves  but  vortex-trapped  near-inertial  internal  waves  (Kunze  et 
al.  1995).  Maximum  velocities  of  the  diurnal  motion  over  Fieberling 
Guyot's  summit  plain  are  ~  15  cm  s“l  some  50-100  m  above  the  bottom  over 
the  summit  plain.  This  signal  decays  •apidly  with  radial  distance  away  the 
seamount.  At  any  instant  and  depth,  the  horizontal  velocity  of  the  diurnal 
motion  forms  a  unidirectional  jet  over  the  summit  which  turns  clockwise  in 
time.  This  horizontal  structure  is  consistent  with  either  the  gravest 
seamount-trapped  mode  (Brink  1990)  or  the  gravest  vortex-trapped  internal 
wave  mode  (Kunze  et  al.  1995).  The  velocity  vector  also  turns 
counterclockwise  with  depth  over  the  entire  200-m  thick  layer  of  intensified 
motion  with  a  vertical  turning  wavelength  Xz  =  250  m.  If  interpreted  as  a 
vortex-trapped  near-inertial  wave,  it  implies  upward  energy  propagation.  It 
contradicts  (i)  Brink's  (1990)  predicted  turning  of  only  ^  90°  confined  to 
within  50  m  of  the  bottom  due  to  bottom  friction  and  a  weak 


counterclockwise-propagating  component,  and  (ii)  the  turning  wavelength 
for  a  time -dependent  nonlinear  benthic  Ekman  layer  (e.g.,  Brink  1995b). 
Beckmann  and  Haidvogel's  (1996)  numerical  simulation  of  Fieberling 
reproduces  most  of  the  above  observed  features  though  their  simulated 
vortex  is  400-  rather  than  200-m  thick. 

Elevated  turbulence  dissipation  rates  e  10-7  w  kg-1  coincide  with  the 
vortex  cap  and  the  diurnal  oscillation,  sufficient  to  drain  the  wave  and  vortex 
in  less  than  three  days  if  they  are  the  energy  sources  for  the  turbulence.  This 
highlights  the  strongly  forced  nature  of  motions  atop  the  seamount's  summit 
plain;  without  the  barotropic  tide  continually  pumping  energy  into  the 
system,  it  would  rapidly  run  down.  The  corresponding  eddy  diffusivities  are 
of  order  10  X  10-4  m2  s-h  These  are  two  orders  of  magnitude  larger  than 
those  deduced  for  the  surrounding  ocean,  suggesting  that,  locally  at  least, 
seamoimts  are  important  sites  for  diapycnal  exchange. 


2.  TRAPPED  Wave  models 


In  a  stratified  fluid,  the  barotropic  diurnal  tide  impinging  on  a 
seamount  will  generate  baroclinic  oscillations  in  order  to  satisfy  the  no¬ 
normal-flow  boundary  condition.  Three  extant  models  may  describe  the 
intensified  baroclinic  diurnal  motions  observed  above  the  summit  plain  of 
Fieberling  Guyot  —  near-resonant  seamount-trapped  topographic  waves 
(Brink  1989),  forced  (off-resonant)  evanescent  topographic  fluctuations 
(Brink  1990;  Eriksen  1995b)  and  vortex-trapped  near-inertial  internal  waves 
(Kunze  et  al.  1995).  We  will  refer  to  slightly  off-resonant  forced  oscillations 
as  waves  here  even  though  they  do  not  exactly  satisfy  the  dispersion  relation 
or  consistency  relations  for  free  waves.  These  three  types  of  waves  share 
many  properties  in  common  and  their  observable  differences  are  subtle.  We 
caution  that  aU  of  these  models  are  approximate.  The  topographic  wave 
models  neglect  background  mean  shear  and  the  vortex-trapped  model  treats 
shear  crudely. 

2.1  Seamount- trapped  topographic  waves 

Brink  (1989)  formulated  a  model  for  free  seamount-trapped  waves  on 
tall  seamounts  in  a  stratified,  rotating  fluid,  extending  the  barotropic  models 
of  Rhines  (1969),  Huthnance  (1974),  Hunkins  (1986)  and  Chapman  (1989). 
These  waves  are  bottom-trapped  topographic  waves  (Rhines  1970)  wrapped 
around  a  seamount  so  that  their  along-isobath  wavenumber  is  quantized. 

For  the  stratification  profile  and  geometry  of  Fiebcrlmg  Guyot,  Brink 
(1990)  showed  that  the  gravest  (highest  frequency)  seamount- trapped  mode 
had  a  slightly  subinertial  frequency  (29-h  period)  and  was  nearly  resonant 
with  the  diurnal  tide.  This  mode  evanesces  rapidly  away  from  the  bottom 
with  a  vertical  decay  scale  of  ~  50  m  and  a  radial  decay  scale  of  a  few 
kilometers  (radial  mode  0).  The  first  azimuthal  mode  is  most  likely  to  be 
excited  by  largescale  forcing  because  its  flow  direction  is  the  same  on  aU 
sides  of  the  seamount  (Brink  1990).  The  direction  of  the  resulting  jet  rotates 
clockwise  in  time  and  propagates  clockwise  around  the  seamount. 


2.2  Forced  evanescent  waves 


Brink  (1990)  extended  his  free-wave  work  to  investigate  the  free¬ 
mode  expansion  of  the  on-  and  off-resonance  forced  response.  Laboratory 
experiments  (Codiga  1993;  Boyer  and  Zhang  1990;  Zhang  and  Boyer  1993) 
and  numerical  simulations  (Haidvogel  et  al.  1993;  Beckmann  and  Haidvogel 
1996)  have  verified  generation  of  baroclinic  subinertial  seamount-trapped 
oscillations  by  barotropic  forcing  in  stratified,  rotating  fluids.  Forced, 
topographicaUy-trapped  motions  can  be  off-resonant  and  thus  do  not  satisfy 
zero  normal  flow  at  the  bottom;  this  is  necessary  for  them  to  compensate  the 
normal  flow  of  the  forcing.  However,  if  they  are  only  slightly  off  resonance, 
evanescent  fluctuations  closely  resemble  free  topographic  waves.  For 
example,  eddy  viscous/diffusive  damping  alters  the  phase  consistency 
relations  for  seamount-trapped  waves  by  only  a  few  degrees  for  the 
wavelengths  of  interest  (Appendix  B).  Observationally,  Brink's  off- 
resonance  predictions  are  consistent  with  the  diurnal  frequency  and 
horizontal  stmcture  of  the  dominant  velocity  fluctuations  above  the  summit 
plain  of  Fieberling  Seamount  (Eriksen  1991;  Noble  et  al.  1994;  Brink 
1995a). 

One  of  the  few  ways  for  distinguishing  between  gravest  free  or  forced 
seamount-trapped  oscillations  and  vortex-trapped  waves  (section  2.3)  is  with 
the  phase  relation  between  radial  velocity  u/  and  vertical  displacement 
signals.  For  (i)  the  cross-slope  standing  modes  discussed  by  Brink,  (ii) 
topographic  waves  propagating  across  the  slope,  and  (iii)  a  slightly  off- 
resonance  forced  bottom-trapped  response,  u/  is  roughly  out  of  phase  with 
and  Ve'.  This  implies  that  inviscid  seamount-trapped  oscillations  have  near¬ 
zero  radial  heat-flux,  <UrT>  si  0,  and  finite  azimuthal  heat-flux,  <vqT>. 

The  horizontal  velocity  vector  for  cross-slope  standing  seamount- 
trapped  modes  undergoes  180®  reversals,  but  no  turning  with  depth.  A 
benthic  Ekman  layer  will  induce  counterclockwise-tuming-with-depth  near 
the  bottom  and  a  velocity  maximum  above  the  bottom.  Brink  (1995b) 
developed  a  model  for  time-dependent  Ekman  layers  in  the  presence  of 
relative  vorticity.  His  expression  for  the  real  part  of  the  wavenumber  (the 
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part  associated  with  phase  propagation  and  turning  of  the  velocity  vector 
with  depth)  is 


where  6  =  Arctan[(4/(coj  ±/ej^)/(C^wo2)],  wq  ~  0(1)  is  the  nondimensional 
vertical  velocity  at  the  top  of  the  Ekman  layer  associated  with  the  Ekman 
downwelling  that  arises  in  this  geometry,  co/  =  ©£  -  nVo/r  ©£  -  =  Ki  + 

C/2  the  intrinsic  frequency  ,feff=f+  C  the  effective  Coriolis  frequency  , 

©£  the  Eulerian  frequency,  V^e  the  azimuthal  velocity  of  the  vortex  and 
C  =  Vo/r  +  dVe/dr  its  relative  vorticity.  Except  near  the  zero  of  ©  -feff,  the 
wavenumber  can  be  approximated  as 


Relation  (1)  will  be  compared  with  the  observations  in  Fig.  14. 

2.3  Vortex-trapped  near-inertial  internal  waves 

Anticyclonic  vorticity  reduces  the  lower  bound  of  the  internal 
waveband  from  the  planetary  Coriolis  frequency /to  an  effective  Coriolis 
frequency /ejjf  (Weller  1982;  Kunze  1985).  As  shown  by  Kunze  et  al.  (1995), 
for  a  vortex  in  solid-body  rotation  and  a  wave  treated  in  cylindrical 
coordinates,  the  effective  Coriolis  frequency  feff=f+  C  (=  0.5/ for  a 
vorticity  of  -  0.5/as  characterizes  the  core  above  Fieberling,  Fig.  3);  half  the 
modulation  arises  from  interaction  with  the  background  shear  dV^/dr,  and 
half  from  advection  (Doppler-shifting)  by  the  background  flow  Fe/r 
(Appendix  B  of  Kunze  1985).  The  effective  Coriolis  frequency  is  well  below 
the  diurnal  frequencies  Ki  =0.933/=  1.87/ej5^and  Oi  =  0.865/=  IJSfeff, 
permitting  diurnal  frequencies  to  propagate  as  near-inertial  internal  waves 
within  the  confines  of  the  negative  vorticity  core.  Like  the  seamount-trapped 
wave,  such  an  oscillation  would  arise  above  Fieberling  to  satisfy  the  no- 


normal-flow  bottom  boundary  condition  (e.g.^  Baines  1982).  Kunze  et  al. 
(1995)  formulated  a  model  for  a  vortex-trapped,  near-inertial  radial  mode 
while  investigating  near-inertial  wave  critical  layers  in  a  warm-core  ring.  In 
this  model,  the  superposition  of  plane  waves  reflecting  off  the  lateral  walls  of 
the  vortex  sets  up  a  radial  mode  which  closely  resembles  a  seamount-trapped 
wave  Li  its  horizontal  stmcture.  Like  the  seamount-trapped  wave,  the  most 
likely  azimuthal  mode  that  will  be  excited  by  largescale  clockwise  rotary 
forcing  is  n  =  -1,  corresponding  to  clockwise  azimuthal  propagation  around 
the  vortex.  The  intrinsic,  or  Lagrangian,  frequency  is  then  o),-  =  co£  -  nVoIr 
Ki  -I-  C/2  =  0.933/-  0.25/=  0.683/  >fejf.  Like  the  gravest  seamount-trapped 
wave,  this  mode  has  maximum  horizontal  velocity  amplitude  near  vortex 
center  and,  at  any  given  instant,  forms  a  nearly  unidirectional  jet  in  the 
vortex  core.  Decay  is  rapid  moving  radially  out  of  the  vortex.  The 
orientation  of  the  jet  rotates  clockwise  in  time. 


Unlike  the  seamount-trapped  topographic  wave,  a  vortex-trapped 
wave  is  an  internal  wave  and  is  therefore  free  to  carry  energy  and 
momentum  upward  into  the  overlying  water  column.  For  a  near-inertial 
wave  of  upward  energy  propagation,  the  horizontal  velocity  vector  turns 
counterclockwise  with  depth.  With  frequency  co/  =  co£  +  C/2  set  by  the 
forcing,  and  lateral  scales  V  dictated  by  the  vortex  dimensions,  the  vertical 
wavelength  is  given  by  the  dispersion  relation 


Xz  -  Xr 


N 


(2) 


where  Xr  is  the  radial  wavelength  (Kunze  et  al.  1995).  Relation  (2)  wiU  be 
compared  to  the  observations  in  Fig.  18.  Since  internal  waves  can  propagate 
freely  only  where  co/  >feff,  that  is,  inside  the  vortex,  they  reflect  off 
horizontal  turning  points  on  the  vortex  sides  to  set  up  radial  modes.  Upward 
propagation  is  blocked  at  a  vertical  critical  layer  at  the  top  of  the  vortex 
where  energy  will  accumulate  until  instabilities  transfer  it  to  dissipative 
turbulence  (Kunze  et  al.  1995).  Thus,  for  a  vortex-trapped  wave,  maxima  in 
energy,  shear  and  turbulent  dissipation  are  expected  above  the  bottom. 


For  upward  and  clockwise  azimuthal  propagation,  the  radial  velocity 
u/  is  180°  out  of  phase  with  the  vertical  displacement  -  TITz,  giving  a 
net  outward  radial  heat-flux  <UrT'>.  In  the  absence  of  irreversible 
processes  (forcing  and  dissipation),  this  wave  flux  will  be  balanced  by  a 
wave-induced  mean  advective  radial  heat-flux.  Including  eddy  diffusivities 
and  viscosities  ranging  up  to  10  X  ICM  m2  s-1,  consistent  with  the 
microstructure  measurements,  does  not  significantly  alter  the  dispersion 
relation  between  frequency,  vertical  wavenumber  and  horizontal 
wavenumber  (2),  nor  the  consistency  relations  between  horizontal  velocity 
components,  u/  and  ve'.  Thus,  it  will  not  alter  the  turning  of  the  velocity 
vector  with  depth.  The  phase  relation  between  a  viscous  wave's  vertical 
displacement  and  the  horizontal  velocities  is  modified  by  at  most  6°  for  the 
vertical  wavelengths  greater  than  100  m  of  interest  here  (Appendix  A). 

2.4  Summary  of  commonalities  and  differences 

In  the  context  of  Fieberling  Guyot,  seamount-  and  vortex-trapped 
waves  share  the  following  properties: 

•  generation  by  the  barotropic  diurnal  tide  to  satisfy  the  zero-normal-flow 
bottom  boundary  condition,  so  their  Eulerian  frequencies 
co£  =  coj  -f-  nVoIr  co/  -  ^2  =  Ki,  Oi. 


•  clockwise  rotation  in  time  of  the  horizontal  velocity  vector  above  the 
summit  plain  such  that  the  radial  and  azimuthal  velocities  are  90°  out  of 
phase. 

•  azimuthal  propagation  clockwise  around  the  seamount  (azimuthal  mode 
n  =  -1). 

•  radial  decay  away  from  the  seamount. 

The  differences  are  more  subtle: 

•  the  vortex-trapped  wave  requires  a  vortex  of  sufficient  strength  to 
permit  a  diurnal  Eulerian  frequency.  The  conditions  for  which  a 


seamount-trapped  wave  can  exist  in  the  presence  of  a  background  shear 
that  makes  part  of  the  domain  effectively  superinertial  are  unclear 
(Codiga  1996). 

•  vortex-trapped  waves  will  carry  energy  and  momentum  up  into  the 
stratified  water  column.  Thus,  they  can  have  energy  maxima  above  the 
bottom  (upward  energy  propagation  implies  counterclockwise  turning 
of  the  horizontal  velocity  vector  with  depth).  Both  free  and  forced 
seamount-trapped  waves  are  bottom-trapped  topographic  waves  so 
might  be  expected  to  have  their  maximum  signal  at  the  bottom  and  little 
turning  of  the  velocity  vector  with  depth.  However,  subinertial  time- 
dependent  benthic  Ekman  layers  also  turn  counterclockwise  with  depth 
and  will  result  in  velocity  maxima  above  the  bottom  (Brink  1995b). 

•  radial  velocity  u/  and  temperature  T  fluctuations  for  vortex-trapped 
waves  propagating  upward  and  clockwise  around  the  vortex  are  in  phase 
to  produce  a  net  outward  radial  heat-flux  <m^T'>  (in  the  absence  of 
forcing  and  dissipative  processes,  the  divergences  of  this  heat-flux  and  a 
nonzero  vertical  momentum-flux  <ve'w'>  drive  a  wave-induced  mean 
circulation,  Ur  and  w,  which  exactly  counterbalances  the  wave  fluxes  so 
there  is  no  net  change  to  the  background  mean  azimuthal  flow  and 
density  fields).  There  is  no  azimuthal  heat-flux  <vqT'>  associated  with 
vortex-trapped  modes.  For  cross-slope  standing  or  propagating  inviscid 
seamount-trapped  waves,  the  radial  heat-flux  is  zero  and  the  azimuthal 
heat-flux  nonzero.  These  phase  relations  are  only  slightly  modified  by 
inclusion  of  eddy  viscosity  and  diffusivity  (Appendices  A  and  B). 


3.  MEASUREMENT  AND  ANALYSIS  TECHNIQUES 
3.1  Data 

Profiles  were  collected  over  the  summit  and  flanks  of  Fieberling 
Guyot  (32°26’N,  127°45'W)  in  March  1991  as  part  of  the  ONR-sponsored 
Topographic  Interactions  program  to  investigate  the  impact  of  seamounts  on 
physical  and  biological  oceanography.  Time-series  and  surveys  were  made 
with  the  High-Resolution  Profiler  (HRP;  Schmitt  et  al  1988),  Sippican 
expendable  current  profilers  (XCP;  Sanford  et  al.  1982;  1993),  and  Sippican 
T-5  (1600  m)  XBTs.  At  the  seamount's  latitude,  the  Coriolis  frequency  is 
7.8  X  10-5  s-1  (22.4  h)  so  both  the  Ki  (24.0  h,  0.933/)  and  Oi  (25.9  h, 

0.865/)  diurnal  tides  are  subinertial.  Barotropic  diurnal  tides  in  this  region 
propagate  northwest  (Schwiderski  1981a,  b). 

Fieberling  Guyot  is  the  northwestmost  member  of  the  Fieberling 
seamount  chain.  For  the  motions  above  the  summit  plain  discussed  here,  it 
appears  to  be  dynamically  isolated  from  its  nearest  neighbor,  Fieberling  H, 
40  km  to  the  southeast  This  is  consistent  with  the  conclusion  of  Roden 
(1991)  that  the  seamount  chain  affects  background  impinging  currents  as  a 
group  only  at  abyssal  depths.  Fieberling  Guyot  rises  from  an  abyssal  plain  of 
4000-4500  m  depth  to  a  summit  plain  at  500-700  m  depth  (Fig.  1).  By 
definition,  a  guyot  has  a  flat  top  (apart  from  a  narrow  pinnacle  southwest  of 
the  geometric  center  of  the  summit  plain  which  attains  440-m  depth)  due  to 
surface  wave  erosion  before  sinking  of  the  peak  below  the  surface.  Radial 
bottom  slopes  on  the  summit  plain  are  ~  0.05  compared  to  flank  slopes  as 
high  as  0.5.  The  abmpt  change  in  slope  near  the  700-m  isobath  will  be 
referred  to  as  the  summit  rim.  The  summit  plain  is  elongated  so  that  the  rim 
radius  is  about  6.5  km  zonaUy  and  half  that  meridionally  (Fig.  1).  Beyond 
the  rim,  Fieberling  Guyot  can  be  described  as  a  Gaussian  of  radial  scale 
12km(Codigal991). 

Other  relevant  measurements  are  two  mesoscale  CTD  surveys 
conducted  in  August  1989  (Roden  1991)  and  April  -  May  1991  (Roden 
1994),  and  year-long  current-meter  mooring  deployments  from  September 
1990  to  September  1991  at  five  sites  on  Fieberling  Guyot  and  two  on  the 
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abyssal  plain  to  the  west  (Wichman  et  al.  1993;  Noble  et  al.  1994;  Brink 
1995a;  Eriksen  1995a,  b).  Moorings  were  located  at  the  center  of  the  summit 
plain  (C,  Fig.  1),  on  the  rim  to  the  northwest  and  southwest  (R2  and  R3),  in 
1500-m  deep  water  on  the  flanks  to  the  northwest  and  southwest  (F2  and 
F3-5)  and  on  the  abyssal  plain  40  km  to  the  northwest  and  southwest.  The  R2 
mooring  was  an  upward-looking  ADCP  while  the  other  moorings  were 
strung  with  discrete-depth  current  meters. 

The  HRP  measurements  include  (i)  day-long  time-series  sampled 
every  3  h  at  the  central  (C)  and  southwest  flank  (F3)  mooring  sites, 

(ii)  repeated  radial  sections  to  the  3000-m  isobath  over  the  western  and 
southeastern  flanks,  and  (iii)  a  survey  around  the  rim  (Fig.  1).  Profiles  on 
the  neighboring  abyssal  plain  20-40  km  from  the  peak  are  described  by  Toole 
et  al.  (1994).  The  stratified  turbulent  benthic  boundary  layer  overlying  the 
flanks  is  examined  in  Toole  et  al.  (1996).  The  freefall  High-Resolution 
Profiler’s  (Schmitt  et  al.  1988)  sensor  suite  includes  a  CTD  and  a  two-axis 
acoustic  current  meter  which  provide  fmescale  vertical  profiles  of 
temperature,  salinity,  pressure  and  horizontal  velocity  relative  to  an 
unknown  but  depth-independent  constant.  Microstructure  sensors  include  a 
high-speed  thermistor,  a  dual-needle  conductivity  probe,  and  airfoil  velocity 
shear  probes  to  estimate  microscale  temperature  and  velocity  dissipation 
rates.  Schmitt  et  al.  (1988)  and  Polzin  (1992)  detail  reduction  of  the  fine- 
and  microstructure  data  to  0.5-m  depth  intervals;  the  microstructure  analysis 
follows  the  methodology  developed  by  Neil  Oakey  at  the  Bedford  Institute  of 
Oceanography.  Dissipation  rate  estimates  assume  three-dimensional  isotropy 
(Gargett  tt  cl.  1984).  The  production/dissipation  balances  of  Osborn  and 
Cox  (1972)  and  Osborn  (1980)  are  used  to  obtain  diapycnal  eddy  diffusivity 
estimates  employing  an  upper-bound  mixing  efficiency  of  0.25  (Oakey 
1982). 

For  the  seamount  cmise,  HRP  was  equipped  with  a  Datasonics,  Inc. 
acoustic  altimeter  to  trigger  release  of  a  descent  weight  when  the  profiler  got 
within  10  m  of  the  bottom.  The  weight-release  criterion  was  effective  over 
the  summit  plain,  but  less  reliable  over  the  steep  flanks  (presumably  because 
of  poor  acoustic  backscatter).  HRP  dives  frequently  terminated  up  to  100  m 
above  the  local  bottom  and  one  impacted  the  bottom. 


Most  of  the  XCPs  were  deployed  along  radial  spokes  from  near  the 
center  of  the  summit  plain  to  the  2500-m  isobath  on  the  flanks  (Fig.  1).  XCP 
profiles  in  water  depths  less  than  1600  m  measured  into  the  bottom.  A 
survey  of  43  profiles  was  conducted  to  the  south  and  east  between  0800  and 
1300  h  on  7  March  (+)  with  19  of  the  stations  re-occupied  between  1330  a.id 
1730  h  (X).  A  survey  of  34  profiles  was  conducted  to  the  north  and  west 
between  0500  and  0900  h  on  8  March  (-I-).  Finally,  a  survey  of  24  profiles 
was  conducted  over  the  summit  plain  a  week  and  a  half  later  between  0500 
and  0830  /i  on  18  March  (o).  Most  of  the  XCP  data  were  collected  within  a 
few  hours  of  the  same  phase  of  the  diurnal  tide.  No  expendable  probes  were 
dropped  within  1  nmi  of  the  moorings  to  prevent  fouling  the  current-meter 
rotors  with  expendable  wire.  As  a  result,  the  southwest  quadrant  was 
undersampled. 

The  expendable  current  profiler  (XCP)  measures  the  horizontal 
velocity  to  within  an  unknown  but  depth-independent  constant  using  the 
voltage  drop  induced  by  seawater’s  motion  through  Earth's  magnetic  field. 
Temperature  is  determined  with  an  XBT  thermistor.  The  maximum 
measurement  depth  of  the  present  model  is  1600  m.  Pre-analysis  processing 
averages  the  data  into  3-m  bins  with  typical  rms  errors  of  ±0.4  cm  S“1  and 
±0.06°C  for  velocity  and  temperature,  respectively  (Sanford  et  al.  1993). 

3.2  Analysis 

Two  features  predicted  to  occur  atop  a  seamount  are  (i)  a  vortex  cap, 
which  should  be  axisymmetric  and  dominantly  in  azimuthal  velocity  ve,  and 
(ii)  a  near-bottom  intensified  azimuthal-mode-one  seamount-  or  vortex- 
trapped  wave  forced  by  the  barotropic  diurnal  tide.  The  trapped  wave 
should  consist  of  a  unidirectional,  slightly  subinertial  jet  over  the  summit 
which  turns  clockwise  in  time  (see  Fig.  16). 

Treating  the  profile  surveys  as  snapshots,  the  vortex  and  wave  signals 
can  be  separated  by  azimuthal  averaging  <*>e.  The  coordinate  system  is 
centered  on  the  seamount  centroid  (32°26.2'N,  127°45.5'W)  as  determined 
by  Codiga  (1991).  The  unknown  barotropic  offset  is  handled  by  arbitrarily 
removing  the  depth-average  over  150-300  m  depth  before  azimuthal 


averaging.  This  emphasizes  perturbations  at  and  below  the  summit  depth. 

To  compensate  for  the  nonuniform  azimuthal  sampling  (Fig.  1),  the 
azimuthal  averages  are  weighted  by  the  sampling  in  the  N,  W,  S  and  E 
quadrants.  This  approach  works  for  velocity  signals  above  the  summit  plain 
but  does  not  handle  the  dipolar  patterns  associated  with  temperature  (vertical 
displacement).  It  also  does  not  work  below  the  summit  plain  depth  where 
topography  constrains  the  flow  to  be  azimuthal  and  dipolar  (Eriksen  1995b). 
While  it  might  be  preferable  to  fit  the  data  to  azimuthal  modes  as  done  by 
Codiga  (1996),  the  present  scheme  captures  the  dominant  signals  above  the 
summit  plain  and  produces  results  indistinguishable  from  those  of  modal 
fitting.  Toole  et  aL  (1996)  explore  averaging  relative  to  local  bathymetry 
orientation. 

Averaging  the  east  and  north  velocities  filters  out  axisymmetric  flow, 
isolating  the  expected  azimuthal-mode-one  diurnal  wave  signal. 
Azimuthally-averaging  azimuthal  and  radial  velocities  filters  out  all 
azimuthal  modes  except  mode  zero,  isolating  axisymmetric  flow.  We  caution 
that  a  vortex  cap  need  not  be  axisymmetric;  numerical  simulations  do  not 
always  find  the  vortex  to  be  axisymmetric  even  in  the  absence  of  an 
impinging  mean  flow  (Haidvogel  et  al.  1993).  Beckmann  and  Haidvogel 
(1996)  find  that  topographic  deviations  from  axisymmetry  and  bottom 
roughness  have  little  impact  on  the  zeroth-order  dynamics  and  the  horizontal 
velocity  signal  but  are  reflected  in  the  mean  density  and  vertical  velocity 
signatures. 


4.  THE  VORTEX  CAP 


4.1  Spatial  structure 

Radial  sections  of  azimuthally-averaged  radial  velocity  <Mr>e 

azimuthal  velocity  <ve>e  are  displayed  in  Fig.  2.  The  profiles  have  been 
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WKB-stretched  to  Nq  -  5.2  X  10-3  s-1  [z  z"  =  |/V(z) -^fz/Vo]  to  remove 

z 

variations  in  the  vertical  scale  and  amplitude  due  solely  to  changes  in  the 
buoyancy  frequency  N.  True  depths  are  indicated  along  the  right  axes. 

Radial  velocities  (Fig.  2a)  are  less  than  3  cm  s-1  over  the  summit  plain 
(r  <  8  km).  The  uniformly  weak  radial  flows  suggest  that  our  azimuthal 
averaging  scheme  is  not  biased.  At  some  depths,  e.g.,  580  sm  (stretched 
meters),  there  is  a  suggestion  of  persistent  in-  or  outflow.  There  is  also  a  hint 
of  a  3  cm  s-1  inflow  below  600-sm  depth  over  the  flanks.  We  hesitate  to 
equate  or  contrast  these  with  the  1-2  cm  s-1  radial  outflow  reported  near  the 
bottom  by  Eriksen  (1991)  and  Brink  (1995a)  because  of  (i)  possible 
problems  with  azimuthal  nonuniformity,  (ii)  the  short  duration  of  our 
sampling,  and  (iii)  the  profile  surveys  were  collected  during  a  time  of 
relatively  weak  radial  flow  based  on  the  R2  ADCP  velocities  (see  Fig.  6). 

The  average  azimuthal  velocity  (Fig.  2b)  displays  a  much  more  clearly 
defined  signal.  A  200-sm  thick  layer  of  anticyclonic  flow  overlies  the 
summit  plain.  Maximum  speeds  of  — 10  cm  s-1  are  found  at  the  rim  (radius 
r  =  7  km).  The  vortex  is  bounded  above  by  a  100-m  thick 
3<ve>0/Oz  =  1.5  X  10-3  s-1  =  0.3V  vertical  shear  layer  and  below  by  a  200-m 

thick  10-3  s-1  =  0.25V  shear  layer. 

Estimates  of  relative  vorticity  ^  =  9<ve>e/0r  +  <VQ>Qlr  (Fig.  3)  are 
noisier  because  they  involve  derivatives.  The  vorticity  has  values  of 
(-  0.50  ±  0.15)/“ over  the  summit  plain  and  vanishes  over  the  seamount 
flanks.  In  the  vortex  core,  the  effective  Coriolis  frequency 
feff=f+  C  =  0.5/=  271/(44.8  h)  <  Ki,  Oi.  Freeland  (1994)  observed  a 
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Figure  2:  Radial  sections  of  azimuthally-averaged  (a)  radial  velocity  <«^>9 
and  (b)  azimuthal  velocity  <V6>9  with  one  standard  error  indicated  by  dark 
shading.  Nonuniform  azimuthal  sampling  (Fig.  1)  is  compensated  for  before 
averaging.  The  depth  coordinate  (left  axis)  has  been  WKB-stretched  with 
respect  to  Nq  =  5.2  X  10-3  s-F  true  depths  are  shown  along  the  right  axis. 
Numbers  along  the  bottom  axis  correspond  to  the  number  of  drops  going  into 
each  average,  darker  and  lighter  silhouettes  to  extremes  in  the  WKB- 
normalized  bathymetry,  and  open  circles  to  current-meter  sites.  The  radial 
velocity  <Ur>e  (a)  is  everywhere  less  than  3  cm  s-1.  The  azimuthal  velocity 
<V0>9  (b)  exhibits  anticyclonic  flow  of  up  to  10  cm  s— 1  over  the  summit  rim  at 
6-8  km  radius  between  600  and  900  sm  (400-  and  800-m  depth). 
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Figure  3:  Radial  sections  of  azimuthally-averaged  relative  vorticity 

<^>e  =  a<ve>e/3r  +  <ve>6/r  with  one  standard  error  indicated  by  dark  shading. 

Relative  vorticities  of  (-0.50  ±  O.lS/are  found  in  a  200-m  thick  layer  over 
the  summit  plain.  Vorticities  vanish  at  radii  r  >  8  km,  indicatmg  that 
<ve>e  0C'r“l.  Weak  vorticity  inside  2-km  radius  may  reflect  measurement 
uncertainty  in  the  position  of  the  vortex  center. 


40-m  thick  vortex  cap  of  similar  strength  in  ADCP  surveys  atop  Cobb 
Seamount  that  was  steady  on  6-7  day  timescales.  However,  during  a  later 
cmise,  Codiga  (1996)  found  vorticities  weaker  than  -0.1/ over  Cobb.  There 
is  a  hint  that  the  vorticity  may  vanish  inside  2-km  radius  (Fig.  3),  but  vortex 
center  position  uncertainties  are  large  enough  that  this  cannot  be  stated 
conclusively.  The  absence  of  a  positive  vorticity  annulus  outside  the 
azimuthal  velocity  maximum  indicates  that  the  azimuthal  velocity  <ve>0  falls 
off  as  r-1  for  radii  7  <  r  <  10  km.  As  discussed  later  in  connection  with  the 
vortex's  potential  vorticity  anomaly,  isopycnals  are  domed  over  the 
seamount  summit  as  would  be  expected  for  an  anticyclonic  vortex.  Vertical 
displacements  are  at  most  4  m  (see  Fig.  7). 

4.2  Comparison  with  current-meter  observations 

Figures  4  and  5  compare  XCP-derived  vorticities  and  azimuthal 
velocities  with  centered  four-day  current-meter  averages  about  the  time  of 
the  7-8  March  XCP  surveys.  The  current-meter  relative  vorticity  is  inferred 
assuming  solid-body  rotation,  ^  =  2<ve>A>  where  <•>,  denotes  the  four-day 
average.  The  current-meter  data  verify  anticyclonic  swirling  between  400- 
and  600-m  depth  (Fig.  4),  albeit  slightly  stronger.  The  current  meters  and 
profilers  both  find  nearly  solid-body  rotation  inside  7-km  radius  (Fig.  5). 
Outside  the  velocity  maximum,  azimuthal  currents  fall  off  more  slowly  in  the 
current-meter  data  than  in  the  profiler  data.  Mean  flows  shallower  than 
400  m  are  weak  and  show  no  reluctance  to  cross  isobaths  (not  shown). 

Differences  between  profiler  and  current-meter  velocities  likely  arise 
because  the  two  measurements  do  not  sample  the  same  part  of  the  flow.  If  the 
vortex  is  not  perfectly  circular,  is  off  center,  or  is  otherwise  not 
axisymmetric,  differences  in  the  profiler  and  current-meter  estimates  would 
result.  Beckmaim  and  Haidvogel's  (1996)  simulations  produce  a  vortex  with 
bands  of  positive  relative  vorticity  localized  between  R2  and  F2,  and  over  the 
NE  flank  (their  Fig.  8);  at  other  azimuths,  little  or  no  positive  vorticity  is 
present  outside  the  velocity  maximum.  Figures  4  and  5  suggest  that 
referencing  the  XCP  velocities  to  150-300  m  produces  offset  uncertainties  of 
at  worst  a  few  cm/s. 


Figure  4:  Comparison  of  the  vertical  structure  of  the  relative  vorticity  at 
r  =  6.7  km  from  the  XCP  surveys  (solid  curve)  with  four-day  smoothed 
current-meter  vorticities  inside  7-km  radius  assuming  solid-body  rotation 
^  =  2<V0>f/r  (-•-).  Standard  deviations  are  shown  for  each  data  type.  XCP 

vorticities  are  -0.35/ to  -0.40/between  500-  and  650-m  depth,  vanishing 
abmptly  shallower  than  400  m  and  more  gradually  deeper  than  900  m. 
Current-meter  vorticities  are  slightly  stronger  at  -0.5/‘to  -0.6/ 


(cm/s) 


Figure  5:  Comparison  of  the  radial  dependence  of  4-day  smoothed  current- 
meter,  <ve>/  (+)  ,  and  azimuthally-averaged  XCP  azimuthal  velocities,  <ve>e 
(♦ ),  in  the  depth  range  440-570  m.  The  dotted  line  is  a  cubic  least-squares  fit 
setting  <V9>  =  0  at  r  =  0  and  r  =  1 1  km.  The  current-meters  indicate  extreme 
azimuthal  velocities  of —14  cm  S“1  at  r  =  6-7,5  km  with  near  solid-body 
rotation  at  smaller  radii,  XCPs  indicate  a  similar  stmcture  but  shghtly 
weaker  velocities  of  ~  -10  cm  s-h 


4.3  Temporal  behavior 

Fortnightly  modulation  of  vortex  strength  is  apparent  in  2-day- 
smoothed  azimuthal  velocities  from  the  rim  moorings  (e.g.,  spanning  500- 
750  sm  at  R2,  Fig.  6).  This  suggests  that  the  vortex  cap,  with  its  dome  of  cold 
water  over  the  seamount  peak,  is  at  least  partly  sustained  against  benthic 
Ekman  pumping  and  turbulent  decay  by  rectification  of  the  Ki  and  Oi  tides. 

The  R2  current-meter  time-series  also  show  that  the  vortex  abruptly 
turns  on  and  off  every  few  months  (Fig.  6).  Similar  behavior  is  seen  at  R3. 
Using  the  R2  and  R3  moorings  to  obtain  time-series  of  the  vortex  center 
position  reveal  no  evidence  that  the  changes  in  the  vortex  strength  were  due 
to  wobbling.  Brink  (1995a)  argues  that  the  vorticity  modulation  is  a  function 
of  the  impinging  zonal  flow  direction.  His  Figure  13  shows  a  tendency  for 
the  vortex  to  be  present  during  eastward  and  suppressed  during  westward 
farfield  flow.  This  is  contrary  to  p-plane  dynamics  predictions.  McCartney 
(1975)  and  Verron  and  LeProvost  (1985)  find  that  westward  flow  produces 
trapped  vortices  while  eastward  flow  generates  a  Rossby  wave  wake  in 
barotropic  p-plane  simulations  of  flow  past  a  smaU-amplitude  seamount.  The 
HRP  and  XCP  measurements  were  taken  at  a  time  of  transition  from  weak 
eastward  to  weak  westward  impinging  flow  (m  ss±1  cm  s-1).  Since  the 
vortex  flow  exceeds  the  impinging  flow,  the  fluid  it  contains  would  be 
isolated  if  the  flow  was  steady. 

4.4  Potential  vorticity  anomaly 

Though  the  HRP  profiles  reveal  a  200-m  thick  positive  density 
anomaly  (<6ae>  =  0.02  gm  cm-3.  Fig.  7a)  associated  with  the  vortex  cap  (site 
C,  r  ~  1  km)  as  compared  to  flank  profiles  (F3,  r  =  7.5  km),  corresponding  to 
an  upward  vertical  displacement  of  4  m,  it  is  too  weak  to  produce  a 
discernible  buoyancy  frequency  anomaly.  The  relative  vorticity  of  the 
vortex  cap  thus  induces  a  negative  potential  vorticity  anomaly  of 
-0.3fN^  (Fig.  7b)  only  slightly  diminished  by  the  weak  positive  buoyancy 

frequency  anomaly.  The  presence  of  a  potential  vorticity  anomaly  is  not 
consistent  with  pure  Taylor-Proudman  dynamics  which  preserve  potential 
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Figure  7:  (a)  Average  HRP  density  anomaly  <5a6>  (dashed),  buoyancy 
frequency  anomaly  <5N>  (solid)  and  vertical  displacement  <^>  (dotted) 
profiles  at  site  C  (r  ~  1  km)  on  the  summit  plain  compared  to  site  F3 
7  5-km  away  on  the  flanks.  The  4-m  doming  of  isopycnals  above  the  summit 
plam  IS  msufficient  to  produce  a  significant  buoyancy  frequency  anomaly 
(b)  Average  profiles  of  relative  vorticity  C// (dotted)  and  potential  vorticity 
anomaly  [<f+  (solid)  atop  the  summit  plain  (r  <  7  km). 

The  vortex  cap  has  a  negative  potential  vorticity  anomaly 


vorticity.  Viscous  damping  or  bottom  friction  would  produce  a  positive 
anomaly  so  cannot  account  for  the  observed  anomaly.  Irreversible  forcing 
or  dissipative  processes  associated  with  tidal  rectification  must  play  a  role  in 
creating  the  negative  potential  vorticity  anomaly. 

If  features  like  this  potential  vorticity  anomaly  are  occasionally  shed 
off  seamounts,  they  might  explain  the  submesoscale  potential  vorticity 
anomalies  found  downstream  of  Ampere  Seamount  (Kunze  and  Sanford 
1993). 

4.5  Cyclo  geos  trophic  balance 

The  steady  azimuthally-averaged  radial  momentum  equation  can  be 
expressed 


where  the  reduced,  or  kinematic,  pressure  P  =  P/po»  P  is  the  dynamic 
pressure  and  the  neglected  radial  acceleration  dUr/dt  is  at  most  5%  of  the 
Coriolis  acceleration.  The  boxed  terms  are  the  cyclogeostrophic  (gradient- 
wind)  balance  expected  for  steady  inviscid  flow  with  curvature  in  the  absence 
of  strong  forcing,  e.g.,  rings  and  Meddies.  Beclanann  and  Haidvogel  (1996) 
fmd  only  15%  deviations  from  balance  in  their  numerical  simulations.  We 
will  refer  to  the  lefthand  side  of  (3)  as  the  effective  Coriolis  acceleration. 

The  unboxed  terms  on  the  right  represent  momentum-flux-divergence 
forcing  b^’  fluctuating  velocities.  Figure  8  compares  the  right-  and  lefthand 
sides  of  a  radially-integrated  cyclogeostrophic  balance 


'  u 


y  2  ^  3 

^ +/Fe] =  5F  =  \hB  dz^g-^  ^hTdz, 


100  m 


(4) 


making  use  of  the  hydrostatic  balance,  dPIdz  =  B,  and  the  local  7,  oe-relation, 
Bcq/BT,  to  infer  pressure  anomalies  6F  from  the  XCP-derived  temperature 
anomalies  57  relative  to  the  outermost  radial  bin  at  Tq  ~  10  km.  The  vertical 
integral  is  taken  downward  from  the  center  of  our  150-300  m  level-of-no- 


Figure  8:  Comparison  of  the  radially-integrated  mean  effective  Coriolis 

2 

acceleration  +f<VQ>e]-dr  (dotted)  with  the  mean  pressure, 

J 

r 

anomaly  5p  =  l<bT>Q  dz  (solid)  as  a  test  of  the  cyclogeo strophic 

(gradient- wind)  balance  (3)  where  r©  ~  10  km  and  z©  =  200  m.  The  quantities 
are  comparable  above  the  summit  plain  for  radii  r  <  4  km.  But  pressure 
anomaly  and  effective  Coriolis  acceleration  have  different  vertical  structures 
for  r  >  4  km. 


motion.  The  effective  Coriolis  acceleration  terms  result  in  a  negative 
anomaly  between  400-  and  600-m  depth  (dotted  curve  in  Fig.  8).  At  radii 
r  <  4  km,  the  pressure  anomaly  5F  in  this  depth  range  (solid  curve)  is  of  the 
same  sign  and  comparable  magnitude,  so  the  vortex  cap  appears  to  be  in 
cyclogeostrophic  balance.  For  larger  radii,  the  magnitudes  of  the  two  sides 
of  (4)  remain  comparable  but  the  vertical  structure  differs.  The  discrepancy 
becomes  even  greater  below  600-m  depth  (not  shown),  where  a  growing 
negative  pressure  anomaly  has  no  balancing  signal  from  the  effective 
Coriolis  acceleration.  This  discrepancy  is  seen  independently  in  both  XCP 
and  HRP  data,  and  is  also  found  in  Beckmann  and  Haidvogel's  (1996) 
numerical  simulation.  The  horizontal  momentum-flux  divergences  on  the 
righthand  side  of  (3)  cannot  account  for  it,  making  at  most  a  10% 
contribution.  The  vertical  momentum-flux  divergence  contribution  might 
be  estimated  from  <Ur^>  ~  2  x  10-3  ^2  s-2  and  the  bottom  slope  Ooby 
recognizing  that  there  can  be  no  flow  into  the  bottom  so  the  ratio  w'lu/  =  a  is 
identical  to  the  bottom  slope  Oq  at  the  bottom,  and  that  this  constraint  relaxes 
away  from  the  bottom.  For  the  sake  of  argument,  assume 
a  =  ao-exp[(z  -  2^,)/(z^  -  zq)]  relaxing  on  a  depth  scale  z^  -  zq  over  which  the 
imbalance  is  observed,  where  zb  is  the  water  depth  on  the  flanks  and  zq  the 
depth  of  the  summit  plain.  Then 


'd<u/w‘> 


dz 


■dr  = 


'd(a<Ur'^>) 


9z 


dr 


Plugging  numbers  into  these  assumptions,  the  vertical  momentum-flux 
divergence  is  only  a  few  percent  of  the  imbalance  so  is  unable  to  account  for 
the  observed  discrepancy  at  r  >  4  km.  Reasonable  alternative  relations  for  a 
yield  similar  results. 

The  sole  remaining  plausible  explanation  for  the  discrepancy  is  that  the 
azimuthal  sampling  is  inadequate  to  filter  out  the  dipolar  horizontal 
temperature  structure  of  the  diurnal  wave.  Over  the  flanks,  there  are  no 
XCP  data  in  the  SW  quadrant.  In  an  extreme  case,  this  would  make  it 
impossible  to  distinguish  between  an  azimuthally-uniform  contribution  and  a 


dipole  with  an  extremum  in  that  quadrant.  We  conclude  that  the  unbalanced 
pressure  anomaly  6F  at  radii  r  >  4  km  and  depths  z  >  400  m  is  likely  due  to 

sampling  bias. 


5.  DIURNAL  Shear 


5.1  Temporal  behavior 

Current-meter  time-series  (Genin  et  al.  1989;  Eriksen  1991;  Noble  et 
al.  x994;  Brink  1995a)  have  established  that  there  is  a  ±15  cm  s~l  diurnal 
oscillation  atop  the  summit  plain  of  Fieberling  Guyot  with  large 
contributions  from  both  Ki  (0.933/)  and  Oi  (0.865/)  tidal  frequencies.  Both 
are  KX)  times  more  energetic  than  the  barotropic  tides  in  the  far  field,  with 
Ki  about  four  times  more  energetic  than  Oi.  Beating  of  these  two  tidal 
constituents  produces  a  subharmonic  fortnightly  cycle  and  a  harmonic  at  M2 
(of  magnitude  BT  -  0.28°C).  The  dominance  of  the  diurnal  frequency  is 
confirmed  by  the  day-long  HRP  profile  time-series  at  C  and  F3  though  the 
tidal  constituents  could  not  be  separated.  Grayscale  HRP  time-series  at  sites 
C  and  F3  (Fig.  9)  reveal  clockwise  rotation  of  the  velocity  vector  in  time  and 
counterclockwise  rotation  with  depth  (downward  phase  propagation) 
between  300-  and  600-m  depth.  Dropped  lag  coherences  for  vertical 
wavelengths  Xz  =  32-128  sm  imply  a  period  of  23.6  h  (0.95/).  The  dominant 
diurnal  signal  in  the  flank  profile  time-series  is  isolated  by  several  hundred 
meters  from  the  bottom. 

Vertical  isopycnal  excursions  of  ±20  m  are  observed  in  the  center  and 
flank  HRP  time-series  (Fig.  9)  with  similar  vertical  and  temporal  scales  as 
those  found  in  velocity.  Their  downward  phase  propagation  over  the  flank 
leads  us  to  conclude  that  the  step-like  stmctures  seen  in  individual 
temperature  and  salinity  profiles  are  produced  by  diurnal  straining,  not 
turbulent  mixing.  At  site  C  on  the  summit  plain,  the  vertical  displacement 
field  is  dominated  by  semidiurnal  harmonics. 

5.2  Spatial  structure 

Azimuthally  averaging  the  east  and  north  velocities  in  the  XCP 
snapshots  filters  out  axisymmetric  flow  associated  with  the  vortex  and 
isolates  azimuthal-mode-one  stmcture  associated  with  the  diurnal  oscillations 
above  the  summit  plain.  DiumaUy  backrotating  the  profiles  to  a  common 
time  before  averaging  did  not  capture  significantly  more  variance  because 
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the  XCP  profiles  were  all  collected  within  a  few  hours  of  the  same  phase  of 
the  diurnal  tidal  forcing. 

Radial  sections  of  azimuthally-averaged,  WKB -normalized  east  and 
north  velocity,  <u>q  and  <v>e,  are  displayed  in  Fig.  10.  Because  of  the  fixed 
vertical  resolution  of  the  measurements,  WKB-stretching  produces  higher 
stretched  vertical  resolution  where  N  is  small,  accounting  for  the  smaller 
apparent  vertical  scales  deep  than  shallow.  A  150-sm  thick  layer  of 
intensified  east  and  north  velocity  (±  15  cm  s-1)  is  evident  immediately  above 
the  summit  plain  (550-700  sm,  350-500  m).  The  signal  decays  on  3-4  km 
radial  scales  away  from  the  summit  but  is  stiU  visible  at  10-km  radius. 

The  efficacy  of  the  averaging  scheme  in  isolating  the  bulk  of  the 
variance  above  the  summit  plain  is  demonstrated  in  Fig.  1 1.  Between  550- 
and  650-sm  depth,  the  variance  in  the  average  (a)  exceeds  that  in  the  residual 
(b)  by  a  factor  of  two.  At  other  depths,  the  residual  contains  more  variance 
than  the  average,  indicating  that  azimuthal-mode-one  structure  dominates 
only  immediately  above  the  summit  plain.  Maximum  average  velocities 
exceeding  15  cm  s-1  lie  between  600  and  650  sm  (400  and  450  m),  about  100 
m  above  the  bottom  inside  2-km  radius.  Velocity  maxima  about  170  m  above 
the  bottom  at  ~  450-m  depth  are  seen  at  the  rim.  Maxima  above  the  bottom 
could  be  consistent  with  either  the  vortex-trapped  wave  scenario  or 
seamount-trapped  oscillations  with  a  benthic  Ekman  layer. 

The  relatively  low  variance  in  both  the  average  and  residual  between 
150-  and  500-sm  depth  is  a  consequence  of  assuming  a  level  of  no  motion  in 
that  depth  range.  This  will  result  in  excess  variance  below  the  summit  depth, 
particularly  in  the  residual  (Fig.  lib). 

The  diurnal  shear  layer  (Fig.  10)  is  present  in  both  east  and  north 
velocities  but  is  out  of  phase  in  the  vertical  in  the  sense  that  the  velocity 
vector  turns  counterclockwise  with  depth.  This  can  be  seen  most  clearly  in 
Fig.  12,  which  displays  a  radial  slice  of  the  orientation  of  the  horizontal 
velocity  vector,  6  =  Arctan(<v>6/<w>0),  in  the  depth  range  500-800  sm 
(300-600  m)  where  the  signal  is  most  energetic.  The  velocity  vector 
orientation  0  increases  monotonicaUy  with  depth  by  240°  in  150  sm.  This 


r  (km) 

0  2  4  8  10  12  (a) 


Figure  10;  Radial  sections  of  azimuthaUy-averaged,  WKB-normalized  east 
velocity  <u>e  (a)  and  north  velocity  <v>0  (b)  over  the  summit  plain  and 
flanks.  Shading  about  the  profiles  indicates  one  standard  error.  WKB- 
stretched  depths  are  indicated  along  the  left  axis  and  true  depths  along  the 
right  axis.  Light  and  dark  silhouettes  indicate  extremes  of  bathymetry,  open 
circles  the  locations  of  current-meters,  and  numbers  at  the  bottom  the 
number  of  drops  going  into  each  average.  A  150-sm  thick  layer  of 
10-15  cm  s-1  velocities  lies  50-100  sm  above  the  summit  plain  and  rim, 
evanescing  radially  away  from  the  seamount.  Extrema  in  <u>q  and  <v>9  are 
not  at  the  same  depths  but  offset  in  the  sense  that  the  velocity  vector  turns 
counterclockwise  with  depth  (see  Fig.  12  and  13).  There  is  also  a 
horizontally-coherent  wave  of  Xz  =  50  sm  in  <v>6  but  not  <u>q  between 
200-  and  400-sm  depth. 


Figure  11:  Radial  sections  of  WKB -normalized  (a)  average  Cartesian  kinetic 
energy  <KEuv>&  —  (<k>q^  +  <v>^)l'l  and  (b)  residual  kinetic  energy 
KE  uv  =  <u'^  +  v'2>e/2.  The  average  kinetic  energy  (a)  has  a  peak  exceeding 
100  cm2  s-2  between  600-  and  700-sm  (400-500  m  depth),  50-100  sm  above 
the  summit  plain.  The  residual  kinetic  energy  (b)  has  more  uniform  overall 
energies  with  no  well-pronounced  peak  above  the  summit,  demonstrating  the 
efficacy  of  the  averaging  scheme  in  isolating  the  feature  over  the  summif 
plain.  The  larger  residual  energy  below  500-sm  depth  is  due  to  removal  of 
150-500  sm  velocities  from  the  profiles  before  azimuthal  averaging. 

Residual  energy  dominates  over  the  flanks  below  the  summit  depth. 


Figure  12:  Radial  section  of  the  orientation  of  the  azimuthally-averaged 
Cartesian  horizontal  velocity,  0  =  Arctan(<v>0/<u>e),  in  the  shear  layer 
(Fig.  10).  The  orientation  0  is  only  contoured  where  the  horizontal  kinetic 
gj^0j-gy  gxceeds  15  cm2  s~2_  The  increase  of  the  orientation  0  by  240  from 
580-^to  720-sm  depth  is  a  consequence  the  velocity  vector  turning 
counterclockwise  with  depth.  The  corresponding  vertical  wavelength  is 
~  210  sm.  The  sense  of  the  turning  is  consistent  with  (i)  near-inertial  waves 
of  upward  energy  propagation  or  (ii)  subinertial  time-dependent  bottom 
Ekman  spirals.  The  orientation  0  also  increases  by  30  in  6-km  radius. 


uniform  vertical  gradient  is  also  seen  in  the  upward-looking  ADCP 
measurements  in  the  bottom  50  m  at  R2  (Fig.  13).  XCP  profiles  show  that  it 
extends  over  a  depth  range  four  times  greater  than  predicted  by  Brink’s 
(1990)  model  (Fig.  13).  While  the  depth  of  the  velocity  reversals  in  Brink's 
seamount-trapped  wave  model  varies  azimuthally  and  in  time,  the  smooth 
turning  evident  in  the  averages  (Figs.  12  and  13)  is  also  typical  of  individual 
XCP  and  HRP  profiles;  no  180°  reversals  were  observed. 

If  interpreted  as  a  near-inertial  internal  wave,  the  orientation  6  is  a 
measure  of  wave  phase.  The  observed  gradient  Oe/Oz  (s  kz)  corresponds  to  a 
vertical  wavelength  =;  210  sm.  The  counterclockwise  turning  with  depth 
implies  downward  phase  and  upward  energy  propagation.  The  orientation 
also  increases  by  30°  in  5-7  km  radius  (Fig.  12),  consistent  with  a  radial 
wavelength  of  ~  80  km  for  planar  geometry.  For  an  internal  wave,  energy 
propagates  along  lines  of  constant  phase  in  the  vertical  plane  so  the 
orientation  contours  would  be  interpreted  as  outward  as  well  as  upward 
energy  propagation. 

The  dependence  of  the  vertical  turning  wavelength  Xz  (=  'lulkz)  on 
relative  vorticity  can  be  inferred  from  the  year-long  R2  ADCP  mooring 
records.  Fits  at  each  depth  bin  of  a  diurnal  oscillation  plus  a  mean  were  made 
to  two-day  segments  of  the  velocity  records  at  depths  of  660-750  sm 
(450-570  m).  The  relative  vorticity  was  estimated  from  the  two-day-average 
mean  azimuthal  velocity  assuming  solid-body  rotation,  C  =  2<V0>,/r.  The 

vertical  turning  wavelength  was  inferred  from  vertical  linear  fits  to  the 

orientation  of  the  diurnal  velocity  signal  weighted  by  the  clockwise-in-time 

N 

diurnal  kinetic  energy,  that  is,  minimizing  X^^cw,(0i  -  -  0o)^  with 

/  =  1 

respect  to  kz  and  Gq.  Examination  of  individual  orientation  profiles  showed 
them  to  be  generally  linear  with  depth  in  the  counterclockwise-with-depth 
sense,  particularly  for  vorticities  stronger  than  -  0.3/. 

The  results  are  summarized  in  Fig.  14  which  displays  the  inferred 
vertical  turning  wavelength  Xz  and  clockwise-in-time  diurnal  kinetic  energy 
versus  relative  vorticity  C/.  The  vortex-trapped  wave  model  prediction  (1) 
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Figure  13:  Comparison  of  the  velocity  orientation  A©  as  a  function  of  WKB- 
stretched  depth  from  Brink's  (1990)  seamount-trapped  wave  model  (dashed), 
the  rim  R2  ADCP  mooring  (solid),  and  XCPs  near  the  center  and  on  the  rim 
of  the  summit  plain  (dotted).  XCP  and  ADCP  gradients  are  roughly 
consistent  in  the  overlapping  depth  range  680-780  sm.  The  model's  bottom 
spiral  is  also  roughly  consistent  with  the  turning  in  this  depth  range. 

However,  the  XCP  profiles  indicate  that  the  phase  gradient  Qz  extends  to  580- 
sm  depth  while  the  seamount-trapped  wave  model  indicates  uniform  phases 
above  700  sm  until  620-sm  depth,  then  a  180°  phase  shift  (velocity  reversal) 
at  shallower  depths. 


Figure  14:  (a)  Clockwise-in-time  diurnal  kinetic  energy  from  2-day  fits  of 
diurnal  sinusoids  to  the  rim  ADCP  records.  The  energies  are  binned  with 
respect  to  relative  vorticity.  Larger  diurnal  signals  and  strong  negative 
vorticities  are  associated  with  each  other,  (b)  Dependence  on  relative 
vorticity  of  vertical  wavelengths  Xz  inferred  from  the  counterclockwise 
turning  of  the  diurnal  horizontal  velocity  vector  with  depth  (small  dots  in 
middle  panel).  Large  dots  are  bin  averages  with  95%  confidence  limits. 

Solid  and  dotted  curves  are  the  vortex-trapped  wave  vertical  wavelength  (2) 
for  fixed  radial  wavelength  Xr  =  30  km,  and  the  larger  root  for  the  time- 
dependent  benthic  Ekman  layer  turning  wavelength  (1)  using  an  eddy 

viscosity  Ve  =  10  X  10-^  m^  s'L  The  average  observed  wavelength  of 
250  ±  20  m  matches  the  vortex-trapped  wave  model  prediction  for  vorticities 
stronger  than  -0.35/.  (c)  The  intrinsic  frequency  co,-  =  Ki  +  (solid)  and 
effective  Coriolis  frequency /ejy=/-i-  ;  (dotted)  as  functions  of  relative 
vorticity. 


for  a  fixed  radial  wavelength  Xr  =  30  km  (from  Fig.  17)  is  displayed  as  the 
solid  curve.  The  larger  root  of  the  time-dependent  benthic  Ekman  layer 
turning  wavelength  (2)  is  shown  as  the  dotted  curve  using  an  eddy  viscosity 
Vg  =  10  X  10-^  m2  s-1  consistent  with  the  microstructure  data  and  a  turbulent 

Prandtl  number  of  one.  The  Ekman  turning  wavelength  has  a  maximum  at 
an  intrinsic  frequency  co,-  =  o}£  +  ^/2  -feff-f  +  C>  or  ^  =  2(coki  -/)•  For 
vorticities  C  <  -  0-35/,  the  observed  bin-averaged  vertical  turning 
wavelengths  Oarge  dots)  are  250  ±  20  sm,  consistent  with  the  vortex-trapped 
wave  model  (solid  curve)  and  an  order  of  magnitude  larger  than  the  time- 
dependent  Ekman  layer  model  (dotted  curve).  To  reproduce  a  vertical 
wavelength  of  250  sm,  the  Ekman  model  would  require  an  eddy  viscosity  of 
1000  X  10-^  m2  s-1,  corresponding  to  a  viscous  decay  time  kz'^IVe  of  only 

4  h.  For  vorticities  of  -0.35/ to  -0.1/,  the  observed  turning  wavelengths 
remain  at  ~  250  sm  while  both  models  predict  smaller  values.  For  vorticities 
of -0.1/ to  0.0,  a  vortex-trapped  wave  cannot  exist  and  the  observed  Xz  bears 

some  resemblance  to  a  time-dependent  Ekman  layer  (dotted  curve). 
However,  we  note  that  there  are  only  a  few  scattered  data  points  and  little 
clockwise-in-time  kinetic  energy  associated  with  the  diurnal  signal  (Fig.  14a) 
at  these  vorticities.  This  reflects  the  tendency  for  weak  diurnal  signals  to 
accompany  weak  vorticities  due  to  fortnightly  modulation  of  both  signals  in 
sync  with  the  fortnightly  beating  of  the  barotropic  Oi  and  Kj  diurnal  tides. 

Vertical  wavenumber  spectra  from  100-900  sm  (50-800  m)  depth 
shown  in  Fig.  15  affirm  the  rotary  character  of  the  signal  with  depth.  The 
shear  layer  over  the  summit  plain  appears  as  a  peak  at  a  vertical  wavelength 
Xz  =  256  sm  in  the  r  =  0-3  km  counterclockwise-with-depth  (CCW) 
spectrum,  most  pronounced  in  the  variance-preserving  spectrum.  There  is 
no  corresponding  enhancement  in  the  clockwise-with-depth  (CW)  spectrum 
at  this  wavelength.  The  CW  spectrum  is  redder  (Jcz~^^)  than  the  GM  model 
(Garrett  and  Munk  1975;  Cairns  and  Williams  1976)  but  has  comparable 
levels  for  Xz  >  50  sm.  At  radii  greater  than  10  km,  there  are  no  significant 
peaks  in  either  the  CCW  or  CW  spectra.  Both  have  the  same  spectral  slope 
but  slightly  lower  (0.7 GM)  levels.  Spectral  slopes  are  steeper  ikz~^)  for 
vertical  wavelengths  Xz  <  50  sm. 
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Figure  15:  Vertical  wavenumber  spectra  of  rotary  kinetic  energy  from 
100-900  sm  (50-800  m  depth)  for  drops  at  0-3  km  (upper  two  panels)  and 
9-12  km  radius  (lower  two  panels).  Spectra  are  presented  in  both  log-log 
(left  panels)  and  variance-preserving  (right  panels)  formats.  The  dotted 
curves  correspond  to  the  GM  model  spectrum.  The  0-3  km  spectra 
emphasize  the  depth-rotary  nature  of  the  shear-layer  signal  above  the  summit 
plain.  They  are  dominated  by  a  counterclockwise-with-depth  (CCW)  peak  at 
Xz  =  256  sm.  There  is  no  corresponding  enhancement  of  CW  energy.  Away 
from  the  seamount,  clockwise  and  counterclockwise  spectra  are  similar  and 
slightly  below  GM  levels. 
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Figure  16  compares  the  horizontal  structure  of  the  amplitude  and 
phase  for  the  Xz  =  256-sm  CCW  component  with  that  for  gravest  seamount- 
or  vortex-trapped  modes.  The  observed  structure  resembles  the  gravest 
mode  in  that  the  largest  vectors  are  found  above  the  summit  plain  and  exhibit 
nearly  uniform  orientation.  Moving  off  the  summit  plain,  the  observed 
vectors  become  weaker  and  of  more  random  orientation.  Comparing  just  the 
radial  decay  of  amplitude,  Fig.  17  displays  an  envelope  of  solutions  for  the 
vortex-trapped  near-inertial  wave  model  constrained  by  the  observed 
buoyancy  frequency  N,  vortex  radius  Tq,  vorticity  C,  and  vertical  wavelength 
Xz  (gray  stippling);  model  amplitude  has  been  normalized  by  the  value  at 
r  =  0.  Further  constraining  the  Eulerian  frequency  to  be  equal  to  the  Ki 
diurnal  frequency  produces  the  black  envelope  of  solutions.  Both  seamount- 
and  vortex-trapped  wave  models,  and  the  XCP  data  (•),  show  the  largest 
signal  inside  the  vortex  and  radial  decay  on  scales  of  a  few  kilometers  outside 
the  vortex.  Since  the  gravest-mode  horizontal  stmcture  for  the  two  wave 
models  is  nearly  identical,  it  cannot  be  used  to  distinguish  which  model  is 
more  appropriate. 

53  Radial  heat-flux 

A  striking  feature  in  the  summit  plain  current-meter  records  is  the 
radial  heat-flux  <Ur'T'>  associated  with  the  diurnal  oscillations  (Brink 
1995a;  Eriksen  1995b).  Johnson  and  Sanford  (1980)  reported  similar 
correlations  between  vertical  displacements  and  outward  velocities  on  the 
flanks  of  Bermuda  which  they  interpreted  as  evidence  of  an  anisotropic 
internal  wave  field.  Radial  velocity/vertical  displacement  correlations  are 
also  seen  in  the  profiler  data  on  the  Fieberling  summit  plain.  Vertical 
wavenumber  coherences  <u/vq>,  <Ur%'>  and  <vq%'>  are  significant  at 
vertical  wavelengths  Xz  =  100-300  sm  for  r  =  6-10  km  (Fig.  18).  At  these 
wavelengths,  u/  and  ve'  are  90°  out  of  phase,  consistent  with  the 
counterclockwise  turning  with  depth  of  the  horizontal  velocity  vector 
(Fig.  12).  Radial  velocity  u/  and  vertical  displacement  -  TITz  are 

roughly  180°  out  of  phase,  implying  an  outward  radial  heat-flux.  It  follows 
that  azimuthal  velocity  ve'  and  vertical  displacement  are  90°  out  of  phase. 
These  properties  are  more  consistent  with  a  vortex-trapped  near-inertial 
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Figure  17:  Radial  structure  of  horizontal  velocity  for  a  vortex-trapped  mode 
(stippled  envelope)  constrained  by  the  observed  range  of  buoyancy 
frequency  N,  core  vorticity  C,  core  radius  Tq  and  vertical  wavelength  Xz-  If 
further  constrained  to  have  the  observed  Ki  diurnal  Eulerian  frequency,  the 
solid  envelope  results,  coi  responding  to  an  inner  radial  wavelength  of 
27  km  and  an  outer  (decay)  radial  wavelength  Vq  =  25  km.  A  seamount- 

trapped  wave  has  identical  radial  structure  (A).  XCP  observations  (f )  Me 
near  the  model  curves. 
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Figure  18;  Cohc’-ence  magnitude  (upper  panels)  and  phase  (lower  panels)  as 
functions  of  vertical  wavenumber  for  r  =  6-10  km.  Phase  is  only  plotted 
where  the  coherence  is  significant  (+)  or  coherence  between  east  and  north 
velocity  components  is  significant  (X).  Phase  relations  for  vortex-  and 
seamount-trapped  waves  are  shown  by  the  light  and  dark  stippling, 
respectively.  Spreading  of  the  phase  relations  by  eddy  diffusivities  or 

viscosities  Kg  and  <  10  X  10“^  m^  s-1  (Appendices  A  and  B)  is  shown  for 
Xz  >  100  m.  Coherence  is  significant  for  vertical  wavelengths  of  100-300  m. 
In  this  band,  radial  and  azimuthal  velocities  are  90°  out  of  phase  (consistent 
with  the  CCW  mining  with  depth),  radial  velocity  u/  and  vertical 
displacement  are  180°  out  of  phase  (corresponding  to  an  outward  radial 
heat-flux),  and  azimuthal  velocity  ve'  and  vertical  displacement  are  90°  out 
of  phase.  These  last  two  phase  relations  are  consistent  with  vortex-trapped 
but  not  seamount-trapped  waves. 
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internal  wave  than  a  seamount-trapped  topographic  wave.  Adding  turbulent 
damping  with  Kp  and  ^  10  X  10-^  m2  s-^,  consistent  with  the  observed 

turbulence  levels,  causes  the  theoretical  phase  relations  to  diverge  from  their 
inviscid  values  (Appendices  A  and  B),  but  not  by  enough  to  explain  the  data 
with  seamount-trapped  fluctuations  (Fig.  12). 

The  radial  structure  of  the  radial  displacement-flux  <UrX>>  shown  in 
Fig.  19  (solid  curve),  resembles  that  of  a  vortex-trapped  wave  (dashed), 
having  a  zero  crossing  near  3-km  radius,  positive  values  inside  that  radius 
and  negative  values  outside.  We  caution  that  Haidvogel  et  al.  (1993)  found  a 
radial  heat-flux  in  numerical  simulations  where  the  forcing  frequency  and 
Coriolis  frequency/ were  too  far  apart  for  a  vortex-trapped  wave  to  exist. 
They  made  no  mention  of  an  azimuthal  heat-flux. 
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Figure  19:  Radial  structure  of  the  real  (co-spectral)  part  of  the  radial 
displacement  flux  UrX  ihe  Xz  =  100-300  m  band.  Dots  and  triangles 
correspond  to  values  from  individual  XCPs  and  HRPs,  respectively.  Solid 
curves  are  radial  bin  averages.  The  dashed  curve  depicts  a  vortex-trapped 
wave  model  prediction. 


6.  FINESCALE  INSTABILITY,  TURBULENCE  AND  MIXING 

Commensurate  with  the  large  amplitude  and  small  vertical  scale  of  the 
diurnal  velocity  signal  above  the  seamount  summit,  the  finescale  gradient 
Richardson  number  Ri  =  N^IVp-  associated  with  these  motions  was 
frequently  ^  0.25  (h’ig.  20).  The  cumulative  probability  distribution  for 
10-m  Ri  above  the  summit  contrasts  markedly  with  that  from  profiles  10  km 
away  (Fig.  21).  In  the  far  field,  only  3%  of  the  10-m  Ri  samples  are  less  than 
0.25,  while  21%  of  the  summit  10-m  Ri  values  lie  below  0.25.  Also,  3%  of 
the  summit  samples  contain  10-m  density  inversions  (10-m  <  0, 

corresponding  to  10-m  Ri<0  in  Fig.  21)  while  no  overturns  are  found  on 
this  scale  in  the  far  field.  Low  Ri  above  the  summit  plain  is  controlled 
principally  by  the  diurnal  shear.  The  shear  of  the  mean  vortex  was  relatively 
small  and  incapable  of  supporting  near-critical  Richardson  numbers  on  its 
own.  However,  its  contribution  to  the  total  shear  is  not  inconsequential;  the 
superposition  of  mean  and  diurnal  shears  modulates  Ri  temporally  and 
azimuthaUy.  Temporal  variability  in  Ri  is  evident  in  the  HRP  time-series  as 
depth-bands  of  low  Ri  propagating  downward  with  the  phase  speed  of  the 
diurnal  velocity  signals.  At  times,  this  results  in  layers  of  Ri  <  0.25  well 
above  the  bottom  as  in  Fig.  20. 

Stability  theory  (Miles  and-Howard  1961)  and  laboratory  studies 
(Thorpe  1978)  suggest  that  flow  at  low  gradient  Richardson  number 
{Ri  <  0.25)  is  susceptible  to  shear  instabilities.  Finescale  density  profiles 
above  the  summit  reveal  frequent  inversions  suggestive  of  the  vertical 
overturning  that  characterizes  this  process  (Fig.  20).  The  histogram  of 
Thoipe  (1977)  overturn  scales  over  the  summit,  where  the  average  buoyancy 
frequency  is  4.5  X  10-3  g-l,  reveals  a  significant  number  of  Thorpe 
overturns  thicker  than  4  m  (Fig.  22).  No  turbulent  patches  with  overturning 
scales  of  this  size  were  found  at  comparable  depths  away  from  the  seamount. 
Nabatov  and  Ozmidov  (1988)  report  similarly  large  overturns  atop  Ampere 
and  Josephine  Seamounts  in  the  North  Atlantic. 

Associated  with  the  low-R/  values  and  thick  vertical  overturns  in  the 
200-m  thick  diurnal  shear  layer  above  the  summit  plain  were  intense 
turbulence  signals  (Fig.  20).  Assuming  3-D  isotropy  as  appropriate  for 
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Figure  20:  A  representative  HRP  profile  from  the  center  of  the  Fieberling 
Guyot  summit  plain.  Panels  include  profiles  of  east  and  north  velocity,  w  and 
V,  with  the  depth-averaged  velocities  removed,  potential  density  ae,  10-m 
linear  least-squares-fit  Richardson  number  Ri  =  N^liuz^  +  Vz^),  Thorpe 
displacements  and  Thorpe  scales  for  each  overturning  patch  (defined  as  the 
interval  over  which  the  sum  of  Thorpe  displacements  is  zeroX  and  turbulent 

kinetic  energy  dissipation  rate  e  averaged  in  0.5-m  bins.  The  elevated 
velocity  signals  below  300-m  depth  are  accompanied  by  more  frequent 
occurrences  of  Ri  <  0.25,  larger  Thorpe  displacements  and  scales,  and  higher 

turbulent  dissipation  rates  e. 
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Figure  21:  The  (a)  probability  and  (b)  cumulative  probability  density 
distributions  of  10-m  Ri  over  the  summit  plain  (thick)  and  10-km  away  from 
the  seamount  (thin).  The  broken  line  demarks  Ri  =  0.25. 
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turbulent  Reynolds  numbers  e/(vA/2)  >  200  (Gargett  et  al.  1984), 
instantaneous  kinetic  energy  dissipation  rates  exceed  e  =  10-6  w  kg-1 
(Fig.  20)  and  temperature  variance  dissipation  rates  exceed  Xj  =  10-*^  s-^ 

within  the  turbulent  patches.  Nabatov  and  Ozmidov  (1988)  similarly 
reported  10-m  thick  patches  of  turbulence  above  and  in  detached  layers  in  the 
vicinity  of  Ampere  and  Josephine  Seamounts  with  turbulent  dissipation  rates 
of  (0.6-40)  X  10-6  Wkg-1. 

Instantaneous  dissipation  rates  of  comparable  magnitude  were  also 
found  in  the  upper  100  m  of  the  water  column  near  Fieberling  Guyot  both  in 
the  surface  mixed-layer  during  periods  of  strong  wind-stress  (e)  and  in  the 
strongly  stratified  pycnocline  immediately  below  it  (e  and  Ij)  (Fig-  20).  The 
signals  just  below  the  base  of  the  mixed  layer  are  associated  with  the  upper- 
ocean  internal  wave  field  and  exhibit  no  systematic  spatial  variation  with 
respect  to  the  seamount.  In  contrast,  between  300-  and  700-m  depth,  the 
turbulent  dissipation  rate  in  the  150-m  thick  diurnal  shear  layer  above  the 
summit  plain  is  consistently  10-100  times  larger  than  that  at  similar  depths 
10  km  away  (Fig.  23a). 

The  verticaUy-averaged  kinetic  energy  dissipation  rate  below  300-m 
depth  above  the  summit  plain  (based  on  profiles  at  or  inside  the  summit  rim) 
is  (3.3  ±  0.1)  X  10-8  w  kg-l.  One  measure  of  the  significance  of  this 
dissipation  rate  is  the  time  it  would  take  it  to  quench  energies  of  the  diurnal 
oscillation  and  vortex  cap.  These  have  average  total  energy  densities  of 
~  5  X  10-3  ni2  s-2  and  2.5  x  10-3  m^  s-2,  respectively,  comprised  mostly  of 
kinetic  energy.  Assuming  these  flows  are  the  energy  source  for  the 
turbulence,  the  observed  average  dissipation  rate  would  drain  the  diurnal 
oscillation  and  vortex  in  about  three  days.  By  comparison,  it  takes  open- 
ocean  turbulence  100  days  to  deplete  a  canonical  Garrett-Munk  internal  wave 
field  (Garrett  1979).  Barotropic  tides  must  be  supplying  energy  at  a  rate 
comparable  to  the  dissipation  rate  to  maintain  the  observed  signals.  This 
further  explains  the  spring-neap  fluctuations  in  both  the  strength  of  the 
vortex  cap  (Fig.  6)  and  diurnal  shear;  this  strongly  dissipative  system  is 
incapable  of  ringing  for  any  appreciable  time.  Irreversible  turbulent 
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Figure  23:  Average  kinetic  energy  dissipation  rates  e  (a)  and  diapycnal  eddy 
diffusivities  Kg  (b)  as  a  function  of  depth  and  radial  distance  from  the 
seamount.  The  0.5-m  e  estimates  are  averaged  into  100-m  thick  vertical  bins 
and  lateral  bins  grouped  (i)  above  the  summit  plain,  (ii)  at  the  summit  rim, 
(iii)  about  the  1500-m  isobath,  (iv)  about  the  2500-m  isobath,  and  (v)  more 
than  10  km  from  the  seamount.  Lateral  spacing  of  the  average  profiles  is  not 
to  scale  and  the  topography  is  schematic.  The  bracketting  95%  uncertainties 
were  computed  with  a  bootstrap  technique  (Efron  and  Gong  1983;  Polzin 
1992)  and  do  not  differ  substantially  from  those  computed  assuming  a  log¬ 
normal  distribution  (Baker  and  Gibson  1987).  Eddy  diffusivities  (b)  were 
inferred  from  the  kinetic  energy  dissipation  rates  following  Osborn  (1980) 
using  a  mixing  efficiency  of  0.25  (Oakey  1982).  Error  bounds  in  (b)  are 
95%  confidence  limits  from  (a). 


processes  clearly  play  a  major  role  in  the  flow  dynamics  above  Fieberling 
Guyot. 

Mirroring  the  radial  patterns  in  dissipation  rate  e,  the  average 
diapycnal  eddy  diffusivities  inferred  from  the  microstmcture  data  increase 
by  two  orders  of  magnitude  approaching  the  seamount  summit  (Fig.  23b). 
Background  diffusivities  of  ~  0.1  X  10~^  m^  s~l,  as  typical  of  the  open-ocean 
pycnocline  (Gregg  1989;  Polzin  et  al.  1995),  were  observed  10  km  away 
from  the  guyot.  The  inferred  eddy  diffusivity,  averaged  over  the  volume 
inside  the  summit  rim  between  300-m  depth  and  the  bottom,  is 
(30  ±  20)  X  ICM  m2  s~l,  corresponding  to  a  vertical  buoyancy-flux  of 
0(10-8  ni2  s-3).  Lueck  and  Mudge  (1996)  report  a  comparable  average 
diffusivity  above  Cobb  Seamount  at  47°N,  170  km  north  of  Fieberling 
Guyot.  Excluding  depths  shallower  than  300  m,  where  thermohaline 
interleaving  is  evident  and  might  enhance  Xj  estimates,  roughly  consistent 
diffusivities  are  estimated  from  temperature  variance  dissipation  rates 
[Kj  =  X7-/(27’z2),  Osborn  and  Cox  1972]  and  kinetic  energy  dissipation  rates 
{Kp  <  y  e/i^,  Osbom  1980)  using  a  mixing  efficiency  y  =  0.25  (Oakey  1982). 

The  layer  of  intensified  turbulence  above  the  summit  plain  extends 
well  into  the  stratified  water  coluinn  so  the  concern  that  near-bottom 
turbulence  might  contribute  little  mixing  because  it  acts  on  already  well- 
mixed  water  (Garrett  1990)  is  not  germaine  here.  For  reference,  the  height 
scale  for  a  steady  Ekman  layer  is  only  8  m  for  an  effective  eddy  viscosity  of 
50  X  10~4  m2  s-1  (assuming  a  turbulent  Prandtl  number  of  one).  At  greater 
depths,  Toole  et  al.  (1996)  and  Eriksen  (1995a)  infer  eddy  diffusivities  of 
(2.0-6.0)  X  1(H  m2  s-1  in  the  stratified  boundary  layer  above  the  flanks  of 

Fieberling.  There,  the  source  of  energy  for  the  turbulent  mixing  appears  to 
be  high-frequency  internal  waves  critically-reflected  from  the  sloping  sides 
of  the  seamount. 
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7.  Summary  and  Discussion 

Fine-  and  microstructure  profiles  collected  over  Fieberling  Seamount 
reveal  a  200-m  thick  layer  of  intensified  activity  overlying  Ae  summit  plain. 
In  this  layer  co-exists  (i)  an  anticyclonic  vortex  cap  of  relative  vorticity 
-0.5/ (Fig.  2b),  (ii)  0(N)  subinertial  diurnal  shear  (Figs.  9  and  10),  and 
(iii)  elevated  turbulence  with  average  kinetic  energy  dissipation  rates 
e  =  (3.3  ±  0.1)  X  1(H  W  kg-1  (Fig.  23a)  and  microstructure-inferred  eddy 
diffusivities  of  Kp  =  (30  ±  20)  X  ICH  m2  s-i  (Fig.  23b).  All  these  motions 
appear  to  be  driven,  at  least  in  part,  by  the  barotropic  Ki  and  Oi  diurnal 
tides.  The  turbulent  dissipation  rate  implies  a  decay  time  of  about  three  days 
for  the  diurnal  oscillation  and  vortex,  so  barotropic  tides  must  be  continually 
pumping  energy  into  the  summit  motions  to  maintain  them. 

7.1  Vortex 

The  vortex  was  roughly  in  cyclogeostrophic  balance  inside  r  =  4  km 
(Fig.  8).  At  larger  radii  and  depths  greater  than  600  m,  the  observed  radial 
pressure  gradient  and  effective  Coriolis  acceleration  do  not  balance.  This 
imbalance  cannot  be  explained  by  contributions  from  momentum-flux 
divergences  and  appears  to  be  an  artifact  of  inadequate  azimuthal  sampling. 
Beckmann  and  Haidvogel  (1996)  find  the  vortex  cap  in  their  simulations  to  be 
in  near  balance. 

The  vortex  cap  appears  to  be  maintained  by  both  (i)  squashing  of 
impinging  background  flow  and  (ii)  tidal  rectification.  The  vortex  strength 
changes  dramatically  every  few  months  (Fig.  6)  depending  on  the  direction 
of  the  impinging  background  flow  (Brink  1995),  suggestive  of  Taylor- 
Proudman  effects.  The  impinging  flow  would  have  to  renew  water  over  the 
summit  plain  every  three  days  to  maintain  a  Taylor  cap  against  turbulent 
dissipation.  This  requires  impinging  flows  strong  enough  to  cross  the 
summit  in  less  than  3  days  time,  or  flows  of  at  least  3  cm  s~l  while  impinging 
flows  are  generally  about  2  cm  S“1  (Brink  1995). 
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Evidence  for  vortex  maintenance  by  tidal  rectification  includes  (i)  the 
short  turbulent  decay  time,  (ii)  vortex  strength  exhibiting  a  fortnightly  cycle 
of  ±0.3/  (Fig.  6),  and  (iii)  the  vortex  core  containing  a  negative  potential 
vorticity  anomaly  (Fig.  7).  Recent  work  by  Codiga  (1996)  indicates  that 
forcing  at  the  bottom  and  dissipation  within  the  water  column  leads  to 
rectification  in  the  right  sense  to  produce  an  anticyclonic  vortex  for 
seamount-trapped  waves.  Beckmann  and  Haidvogel  (1996)  suggest  that 
vortex-trapped  waves  will  lift  an  existing  vortex  up  off  the  bottom.  The 
strong  role  of  rectification  at  Fieberling  raises  doubts  about  past 
interpretations  of  observed  isopycnal  doming  over  seamounts  in  terms  of 
Taylor-Proudman  dynamics. 

Brink  (1995)  argued  that  the  outward  radial  heat-flux  could  be 
maintaining  the  vortex  cap's  dome  of  cold  water  over  the  summit  against 

benthic  Ekman  downwelling  wT z  [w  =  (-  0.08  ±  0.05)  cm  s-1  as  he  inferred 
from  observed  near-bottom  radial  outflow  of  Mr  =  1-2  cm  s-1, 

Tz=  10~2  °C  m-1].  Current-meter  radial  heat-fluxes  are 
<UrT'>  =  0.4°C  cm  s-1  (Brink  1995;  Eriksen  1995b),  an  order  of  magnitude 
higher  than  the  fluxes  found  here  (Fig.  19)  though  the  profiler  surveys  took 
place  during  a  neap  and  do  not  include  barotropic  velocities. 

There  are  subtleties  in  the  rectification  problem  which  frustrate 
inferring  balances  from  incomplete  data.  Consider  steady  {d/dt  -  0), 
azimuthally-averaged  (9/90  =  0)  conservation  of  azimuthal  momentum  and 
buoyancy  within  a  vortex 
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where  ^  =  Vq/v  +  dVQldr  is  the  background  mean  vorticity,  b\  u/,  ve'  and  w' 
are  the  wave  buoyancy,  radial,  azimuthal  and  vertical  velocities,  and  Ur  and  w 
are  the  wave-induced  mean  cross-stream  circulation.  Values  are  deduced 
from  Table  1.  Note  that  the  fortnightly  acceleration  of  the  vortex  azimuthal 
velocity  is  8  Ve/9/  =  5  X  10-7  m  s-2  while  the  fortnightly  rate  of  change  of  the 

vortex  buoyancy  anomaly  is  dBldt  =  4  X  10-12  jq  s-3.  Eriksen's  (1995b) 
estimate  of  <V0'h’'>  assumes  a  balance  between  time  rate  of  change  and 
vertical  advection.  From  (5),  it  seems  to  be  a  factor  of  2-4  too  high  even 
taking  into  account  acceleration  of  the  vortex.  Substituting  (5)  and  (6)  into 
continuity  (7)  and,  recognizing  that  the  constants  of  integration  are  zero  for  a 
spatiaUy-confmed  wave  field,  we  find  that  the  steady  balance  reduces  to 

(8) 

X  10-1  m2  s-1 

where  the  imbalance  can  be  accounted  for  by  the  factor  of  2-4  overestimate 
of  <V0'w'>.  This  relation  holds  exactly  for  inviscid  unforced  vortex-trapped 
waves.  That  is,  wave  fluxes  exactly  balance  wave-induced  mean  advective 
fluxes  to  maintain  a  steadystate;  the  background  field  remains  unchanged  by 
the  waves'  flrx-divergences  unless  irreversible  forcing  or  dissipation  break 
interrelationship  (8).  This  is  a  general  result  for  waves  propagating  in  a 
background  mean  shear  known  as  the  nonacceleration  theorem  (Andrews  and 
McIntyre  1976;  Dunkerton  1980).  Thus,  with  only  partial  measurements  of 
the  flux  divergences,  the  net  impact  of  the  forcing,  fluctuations  and 
dissipation  on  the  background  cannot  be  evaluated. 

Hopes  of  quantifying  wave-induced  mean  accelerations  from 
numerical  models  where  dissipative  processes  are  prescribed  using 
unverified  parameterizations  are  likewise  fmstrated.  Beckmann  and 
Haidvogel's  (1996)  Fieberling  simulation  is  relatively  insensitive  to  the 
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Quantities  for  the  azimuthally-averaged  equations  (5-7). 


variables 

value 

source 

f 

8  X  10-5  s-1 

/+c 

4x  10-5  s-l 

fjf2 

2xia-5s-2 

fz 

10-2°Cm-l 

50  m 

(Fig.  14) 

L 

5  km 

4m 

(Fig.  7) 

B  =  mo 

-8  X  10-5 m  s-2 

<w’b'> 

3  X  10-3  m2  s-l 

(Fig.  23) 

Ur 

1.5x10-2  ms-1 

(Brink  1995) 

w 

-  8  X  10-^  m  s-l 

(Brink  1995) 

<u/Vq'> 

0 

(Fig.  18) 

<VqV> 

0 

(Fig.  18) 

<Ur'b’> 

10-2  in2  s-3 

(Brink  1995;  Eriksen  1995b) 

0.1  X  10-7  m2  s-2 

(Fig.  19) 

<V0'w'> 

10-4  ni2  s-2 

(Fig.  7  of  Eriksen  1995b) 

magnitude  of  bottom  friction  and  epineutral  diffusivities,  but  comparable 
isosigma  mixing  produced  damping  times  of  300  seconds.  They  find  a 

balance  between  UfdBldr  and  in  contrast  to  the  numbers  in  (6),  and 

radial  heat-fluxes  in  the  opposite  sense  to  the  observations,  that  is,  acting  to 
remove  the  cold  dome  (baroclinic  instability).  Without  knowing  the  detailed 
nature  and  spatial  distribution  of  the  forcing  and  dissipative  processes  in 
models,  or  measuring  all  the  flux-divergences  in  the  field,  the  net  rate  of 
change  of  background  dynamic  quantities  cannot  be  determined 
unambiguously. 

7.2  Diurnal  shear 

The  diurnal  shear  above  the  summit  plain  also  displays  a  fortnightly 
modulation  of  Ki  beating  against  Oi  (Noble  et  al.  1994;  Brink  1995).  Its 
jet-like  horizontal  stmcture  rotates  clockwise  with  time  and 
counterclockwise  with  depth  (kz  =  250  sm).  Its  slightly  subinertial 
frequency  (Fig.  9)  and  azimuthal-mode-one  horizontal  stmcture  (Figs.  16 
and  17)  closely  resemble  those  of  a  seamount-trapped  topographic  wave 
(Brink  1989;  1990).  However,  vortex-trapped  near-inertial  internal  waves 
share  these  properties  (Kunze  et  al.  1995).  The  (i)  continuous 
counterclockwise  turning  of  the  velocity  vector  with  depth  (Figs.  12  and  13), 
(ii)  nonzero  radial  heat-flux,  and  (iii)  near-zero  azimuthal  heat-flux  (Figs.  18 
and  19)  more  closely  resemble  a  vortex-trapped  internal  wave  than  either 
free  or  forced  seamount-trapped  topographic  waves.  The  heat-fluxes  are  the 
most  convincing  argument  in  favor  of  a  vortex -trapped  wave  interpretation 
as  subinertial  benthic  Ekman  layers  also  exhibit  counterclockwise  turning 
with  depth  albeit  with  smaller  scales  (Fig.  14). 

We  conclude  that  the  layer  of  intense  diurnal  shear  observed  in  the 
HRP  and  XCP  measurements  atop  the  summit  plain  is  best-described  as  a 
vortex-trapped  near-inertial  internal  wave,  its  subinertial  frequency  allowed 
by  the  strong  negative  vorticity  of  the  vortex  cap.  With  its  frequency 
determined  by  diurnal  tidal  forcing  and  its  horizontal  scale  controlled  by  the 
vortex  size,  the  vertical  wavelength  of  a  vortex-trapped  wave  is  set  by  the 
wave  dispersion  relation  (2). 


The  inviscid  vortex-trapped  wave  model  (Kunze  et  al.  1995)  predicts 
vertical  and  horizontal  structure  for  tidally-generated  diurnal  motions 
consistent  with  the  observations,  but  it  does  not  predict  amplitude.  The 
Richardson  number  and  microstructure  observations  suggest  that  turbulent 
dissipation  driven  by  shear  instability  limit  the  magnitude  of  the  diurnal 
oscillations.  Energy  supplied  by  the  barotropic  tide  impinging  on  the 
seamount  will  raise  the  amplitude  of  the  vortex-trapped  wave  until  unstable 
Froude  number  \Vz\/N  is  achieved.  Subsequent  energy  will  be  dissipated  by 
shear-driven  turbulence  to  sustain  a  sufficiently  unstable  \Vz\/N  to  account 
for  the  dissipation  (Kunze  et  al.  1990).  Thus,  equilibrium  vortex  and 
diurnal  amplitudes  will  be  established  such  that  the  energy  being  supplied  by 
the  barotropic  tide  is  lost  to  turbulent  dissipation  at  the  same  rate.  This 
equilibrium  will  be  set  up  on  a  timescale  of  a  few  days  since  the  dissipation 
timescale  is  three  days;  vortex  and  diurnal  strengths  track  the  spring-neap 
cycle  of  the  barotropic  tide  with  no  discernible  lag.  For  Fieberling  Guyot’s 

N  =  4.5x  10“3  rad  s-1  and  a  vertical  wavelength  Xz  =  210  m,  a  Froude 
number  \Vz\/N  =  2.0  corresponds  to  velocities  of  30  cm  s-1  which  are 
achieved  when  the  diurnal  and  vortex  shear  are  additive. 

How  then  do  we  account  for  the  evanescent  diurnal  signal  found  over 
the  flanks  that  is  coherent  with  motions  above  the  summit  plain  (Eriksen 
1995b).  Vortex-trapped  waves  cannot  exist  on  the  flanks  so  these  must  be 
bottom-trapped  motions.  This  is  consistent  with  the  flank  azimuthal  heat- 
flux  dominating  over  the  radial  heat-flux  (Eriksen’s  Fig.  7).  He  models  them 
as  forced  evanescent  oscillations,  that  is,  off-resonant  forced  topographic 
waves.  The  downslope  phase  propagation  he  observed  then  implies  upslope 
energy  propagation.  This  upslope  energy-flux  must  be  accumulating  and 
dissipating  on  the  summit.  We  speculate  that  it  feeds  energy  into  the  summit 
vortex-trapped  motions. 

Vortex-trapped  internal  waves  forced  by  diurnal  tides  should  be  found 
only  just  poleward  of  30°.  Equatorward  of  30°,  diurnal  frequencies  can 
propagate  as  free  internal  waves  so  will  not  be  trapped.  Poleward  of  43°,  a 
vortex  cap  could  not  be  strong  enough  to  allow  a  finite-scale  free  diurnal 


wave  within  its  confines  (Fig,  24).  This  is  consistent  with  Codiga's  (1996) 
interpretation  of  the  diurnal  signal  in  repeated  ADCP  surveys  over  Cobb 
Seamount  (46°45'N)  as  a  seamount-trapped  wave.  Thus,  while  intensified 
baroclinic  tidal  oscillations  appear  to  be  a  common  feature  of  seamount 
summits,  their  dynamics  appears  to  be  a  strong  function  of  latitude. 

Vorticity-trapped  diurnal  waves  are  also  likely  to  arise  along  the 
continental  shelf  break.  Tidal  rectification  will  drive  poleward  flow  on  the 
eastern  boundary  and  equatorward  flow  on  the  western  boundary.  Negative 
vorticity  associated  with  these  currents  will  trap  subinertial  internal  waves 
driven  by  diurnal  tides  for  latitudes  ~30°-40°,  by  semidiurnal  tides  for 
latitudes  greater  than  70°  and  by  annual  tides  near  the  equator.  The  cross¬ 
slope  heat-flux  <u'T'>  and  vertical  momentum-flux  <v'w'>  associated  with 
these  trapped  oscillations  may  redistribute  or  drive  stronger  along-slope 
currents,  such  as  the  California  Undercurrent  (Hickey  1979). 

7.3  Turbulent  mixing 

Elevated  turbulence  levels  of  e  ^  10-7  W  kg-1  are  found  coincident 
with  the  vortex  cap  and  diurnal  oscillation,  corresponding  to  eddy 
diffusivities  of  0(10  X  1(M  nfi  s-l).  Geometric  scaling  arguments  provide  a 
means  of  assessing  the  local  and  global  importance  of  seamount  mixing 
(Lueck  and  Mudge  1996).  Applying  their  approach  to  Fieberling,  the  layer 
of  intensified  turbulence  atop  the  seamount  has  a  horizontal  area  of 
A5  ~  100  km2.  Therefore,  the  average  inferred  diffusivity  on  the  summit 
implies  a  diapycnal  buoyancy-flux  of 

Ft,  =  KpAsN^  ~  (25  X  10*^  m^  s~l)  X  (10^  m^)  x  (10-5  s-  7)  =  2.5  m^  s-3. 
Substantial  evidence  has  accumulated  in  recent  years  to  support  an  interior 
(i.e.,  away  from  boundaries)  diapycnal  diffusivity  -  0.1  X  ICH  m2  s-1 
(Mourn  and  Osborn  1986;  Gregg  1989;  LedweU  et  al.  1993;  Toole  et  al. 

1994;  Polzin  et  al.  1995).  Thus,  the  interior  ocean  area  required  to  support 
the  same  diapycnal  buoyancy-flux  found  over  the  seamount  summit  is 
2.5xl04km2,  or  250  times  the  area  of  the  summit  plain  mixing  zone. 

Applying  the  same  arguments  on  the  basin  scale,  consider  a  density 
surface  that  intersects  the  intense  mixing  layer  atop  the  seamount.  Ignoring 
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Figure  24:  (a)  Intersects  of  the  tidal  constituent  frequencies  (Neumann  and 
Pierson,  1966)  and  the  Coriolis  frequency/ as  a  function  of  latitude.  Decadal 
through  fortnightly  tides  cluster  within  2°  of  the  equator,  diurnal  tides^ 
straddle  30°  and  semidiurnal  tides  resonate  at  latitudes  greater  than  68°. 

(b)  The  vertical  wavelength  Xz  of  a  Ki  Eulerian  frequency  internal  wave  with 
horizontal  wavelength  30  km.  The  dotted  curve  corresponds  to  a  quiescent 
ocean  where  free  diurnal-frequency  internal  waves  are  not  possible  poleward 
of  30°.  The  solid  curve  corresponds  to  a  vortex  of  core  relative  vorticity 
—0.5/ and  azimuthal  mode  n  —  —1  for  which  diurnal  internal  waves  can  exist 

equatorward  of  42°. 


that  this  surface  may  outcrop  farther  north,  its  surface  area  north  of  20°N  in 
the  Pacific  is  about  3.8  X  10^  km2.  Hypsometric  data  reveals  that  0.014  of 
this  surface  lies  within  300  m  (the  thickness  of  the  intensely  turbulent  layer 
above  Fieberling's  summit  plain)  of  the  ocean  bottom  (see  Fig.  7  of  Toole  et 

al.  1996).  For  summit  values,  we  deduce  a  diffusive  diapycnal 

buoyancy-flux  for  the  interior  N.  Pacific  pycnocline  of 

(3.8  X  1013  m2)  X  (0.1  X  10-4  m2  s-1)  X  (10-5  s-2)  =  4  X  103  m4  s-3.  If  we 

assume  that  the  average  diffusivity  in  the  majority  of  the  near-boundary 
regions  is  comparable  to  that  observed  above  the  summit  plain,  the 
corresponding  near-boundary  buoyancy-flux  is  (0.014)  X  (3.8  X  1013  m2)  x 
(25  X  ICM  m2  s-1)  X  (10-5  s-2)  =  13  X  103  m4  s-3^  three  times  larger  than 
the  interior  buoyancy-flux.  The  total  (interior  plus  near-boundary)  inferred 
buoyancy-flux  is  about  half  what  would  be  obtained  with  a  spatiaUy-uniform 
diapycnal  diffusivity  of  1(H  m2  s-1  acting  over  the  whole  surface.  However, 
only  a  small  fraction  of  the  near-boundary  region  is  occupied  by  seamount 
summits.  Continental  slopes  and  flanks  of  midocean  ridges  and  seamounts 
account  for  most  of  the  near-boundary  area.  Toole  et  al.  (1996)  and  Eriksen 
(1995a)  estimate  effective  diffusivities  of  (2-6)  X  KM  m2  s-1  for  the 

stratified  turbulent  boundary  layer  adjacent  to  the  flanks  of  Fieberling. 
Taking  this  diffusivity  as  more  representative  of  the  near-bottom  average 
reduces  the  near-boundary  buoyancy-flux  to  (1-3)  X  103  m4  s-3,  comparable 
to  the  interior  flux.  Under  this  assumption,  the  interior  plus  near-boundary 
buoyancy-flux  falls  to  ~  6  X  103  m4  s-3,  only  15%  of  what  would  result  from 
a  SpatiaUy-uniform  diapycnal  diffusivity  of  KM  m2  s-1. 

From  the  above  arguments,  we  conclude  that,  on  a  local  scale, 
seamounts  like  Fieberling  Guyot  constitute  very  significant  sites  for 
diapycnal  exchange  as  also  found  by  Lueck  and  Mudge  (1996).  However,  on 
basin  scales,  there  does  not  appear  to  be  sufficient  ocean  area  in  close 
proximity  to  seamount  summits  to  support  average  diapycnal  transports 
equivalent  to  those  that  would  be  associated  with  the  canonical  (Munk  1966) 
KM  m2  s-1  diffusivity  in  the  main  pycnocline.  Ideal  thermocline  theory 
(e.g.,  Luyten  et  al.  1983)  demonstrates  that  diapycnal  transports  of  this 
magnitude  are  not  necessary  in  the  main  pycnocline. 


In  contrast,  direct  mass  and  heat  budgets  for  abyssal  basins  (e.g.,  Hogg 
et  al  1982;  Roemmich  et  al  1995)  do  require  average  diapycnal  diffusivities 
of  0(1(H  m2  s-1).  As  Toole  et  al.  (1996)  discuss,  the  characteristic 
hypsometry  of  deep  ocean  basins  results  in  abyssal  density  surfaces  having  a 
large  fraction  of  their  area  in  close  proximity  to  the  bottom  (their  Fig.  7). 
The  average  diffusivity  on  a  given  density  surface  is  obtained  oy  weighting 
the  high  diffusivities  near  topography  by  the  fraction  of  the  ocean  basin 
proximal  to  bathymetry  on  that  surface. 


^  Ks-As  Ki(Atot-As)  KsAs 

^  ^tot  ^  ^tot  ^  ^tot 


(9) 


Taking  near-boundary  diffusivities  Ks  comparable  to  those  found  atop  the 
summit  plain,  waters  below  3000-m  depth  wiU  have  average  diapycnal 
diffusivity  raised  above  10^  m2  s-1  while,  for  near-boundary  diffusivities 
close  to  observed  flank  values,  depths  below  4500  m  will  have  average 
diffusivities  exceeding  10~^  m2  s~l  due  to  near-boundary  mixing.  Thus, 
while  topographically-induced  mixing  plays  at  most  a  comparable  role  to 
weak  interior  mixing  in  the  main  pycnocline,  it  could  dominate  near  the 
ocean  bottom  to  drive  0(1(H  m2  s-1)  abyssal  diffusivities. 
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Appendix  A;  Viscous  Near-Inertial  Waves 


Consider  the  equations  of  motion  for  hydrostatic  rectilinear  internal 
waves  of  the  form  V  =  ¥o(r)  exp[/(itxA:  +  kzz  -  co/r)],  where  k  -  (kx,  0,  kz)  is 
the  wavevector  chosen  so  the  horizoniai  direction  of  propagation  lies  along 
the  j:-axis  (ky  =  0)  and  co/  the  intrinsic  frequency  damped  by  eddy  viscosity 
and  eddy  diffusivity  Kg, 


-fv  =  -ikxp  +  VekP-u 


-/co,v  +/w  =  Vekz'^v 


0  =  -ikzp  +  b 


-mib  +  N^w  =  K^kPb 


kxU  +  kzW  =  0 


-ib 

t: 


(Al) 


SO  that 


[aii-iVekP]-u-Jv:= 


N^kx^u 

(co,  -  iKekz^)kz'2- 


-[/co/  +  Vekz'^\v  ■^fu  =  0 


-ifu 

(co/  -  iVekz^) 


(A2) 


and  the  buoyancy  can  be  espressed  in  terms  of  the  zonal  velocity  u  as 


(A3) 


(£0|  —  iK.gkz^)kz 


Substituting  the  boxed  expression  for  v  in  (A2)  into  the  first  expression  in 
(A2)  and  assuming  weakly-damped  (Vekz^,  «f)  near-inertial 
(N^kp  «pkp)  waves,  the  dispersion  relation  can  be  approximated  by 


CO/  +  i<5  =  /  + 


N^kx^ 

ifkp-  + 


iV  ekz'^ 


(A4) 


The  effect  of  viscosity  appears  as  a  weak  evanescence  <3  in  time  (or  in  the 
direction  of  energy  propagation  Cg).  Within  the  context  of  this 
approximation,  the  eddy  diffusivity  has  no  impact  on  the  dispersion  relation 
because  near-inertial  waves  have  relatively  weak  buoyancy  signals.  More 
general  viscous  internal  wave  behavior  can  be  found  in  LeBlond  (1966).  For 
Ve  =  10  X  10-^  m2  s-1  and  Xz  =  100-300  m,  VekP  =  (0.4-4.0)  X  10-6  s-1  so  is 
much  smaller  than  the  Coriolis  frequency /=  7.8  X  10-6  s-k  Substituting 
dispersion  relation  (A4)  into  (A3)  and  (A2),  we  obtain  the  following 
consistency  relations  between  the  buoyancy  b,  north  velocity  v,  and  east 
velocity  u  . 


bPkx 

■  fkz 


r  bPkp  (Vg  -  Ke)kz^  1 

7  J“ 


V  = 


The  phase  relations  between  uv,  ub  and  vb  are  then 


(A5) 


^uv  ~  ' 


K 


2flkz^  -  N2kx^ 

^ub  =  ArctanL 


IfkzHVe  - 


]  =  Arctan[' 


(V, 
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r-(Ve  -  Ke)/:z2 

<^v6  =  Arctant - - )J, 


so  that  the  two  horizontal  velocity  components  remain  90°  out  of  phase  but 
the  velocity  and  buoyancy  perturbations  have  slightly  jnodified  phase 
relations  if  the  eddy  viscosity  and  diffiisivity  differ  from  zero  and  each  other 

(eddy  Prandtl  number  Vg/Kg  ^  1.0).  For  eddy  viscosities  and  diffusivities 

Vg,  Kg  <  10  X  10“^  m2  s~l  and  Xz  =  ICX)  m,  the  phase  relations  differ  by  at 
most  3°  from  their  inviscid  values  (Fig.  18). 


Appendix  B:  Viscous  Stratified  topographic  Waves 


Following  Rhines  (1970)  and  Appendix  A,  the  equations  of  motion 
rotated  in  the  vertical  plane  to  lie  parallel  to  topography  for  plane  stratified 

bottom-trapped  topographic  waves  with  eddy  viscous  Ve  and  diffusive  Ke 
forcing  on  gentle  slopes  (a  «  1)  are 

-/(co  -  i^th^)u\\  -yV  =  -ik\p  +  ba 


-/(co  -  iVek±^)v  +  fuw  =  -ikyp 

fay  =  -ik^  +  b  (B 1 ) 


-i(co  -  iKek±^)b  =  -A^awii 


=> 


-iN^a-U\\ 
(co  -  iKek±^) 


^||W||  -I-  kyV  —  0 


where  (mh,  v)  are  the  across-  and  along-isobath  velocities  parallel  to  the 
bottom,  the  flow  normal  to  the  bottom  w±  is  assumed  to  vanish  identically, 
and  (/til,  ky,  k±)  are  the  across-isobath,  along-isobath  and  normal 
wavenumbers  in  a  coordinate  system  lying  in  the  plane  of  the  bathymetry.  It 
has  been  assumed  that  kj_ »  k\\,  ky  in  the  viscous  terms,  consistent  with  the 
observations.  Substituting  the  boxed  expressions  for  buoyancy  b  and  along- 
isobath  velocity  v  into  the  first  three  expressions  of  (Bl), 

-i(co  -  iVek±^)  f_  iN^a'2- 

(cD  -  i'Kek±^)k\\ 


]mii  =  -ip 
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Ltt7 + - — J«ii = -ip 

(CD  -  iKeks?)kx 


or,  equating  the  lefthand  sides, 

m2(v2^^2 

(co  -  iVekx^Xa  -  iKek±^) 


(B3) 


(co  -  iVek±^)((o  -  iK£k±^)k\\k±  -  ifkyk±(oi  -  iKek±^)  -  =  0. 


which  can  be  solved  for  co  and  ^lin  terms  of  A^,/,  a,  V^,  IQ,  ky  and  ^h. 
Assuming  that  the  frequency  co  is  altered  only  slightly  so  that  it  can  be 
expressed  as  a  small  perturbation  from  the  inviscid  solution 
CO  =  Nakylkn  +  5cd,  where  5co  «  Nakylkn,  the  first  equation  in  (B3)  implies 


Naky  i(Ve  +  ^e)kx^ 
2 


(B4) 


where  Nakylkn  ^f»  Veki^,  from  before.  Substituting  (B4)  into  the 
second  equation  in  (B3)  gives 

N‘^a^kyk\\k±  ^Naky  -  ^e)k±^\ 

kj.  - 2 - )-W2t^y  =  0  (B5) 

and,  assuiring  a  solution  for  k±  perturbed  about  the  inviscid  solution 


_  Nkn^  _  2 

'^^-Nah-ifkn 


Nak\\  +  ifkn 
N'^a^k\?-  +  P-kn^  ’ 


(B6) 


leads  to  a  viscous  perturbation 


hk\^  = 


N^iVe  -  K,)/://8 

2aky{Nak\\  -  ifknP  ' 


(B7) 


This  could  blow  up  if  ky  =  0,  invalidating  the  small  perturbation 
approximation.  Based  on  the  observational  results  of  a  diurnal  Eulerian 
frequency  co£  =  Ki  and  an  imaginary  vertical  wavenumber  at  most 
comparable  to  the  real  vertical  wavenumber,  kz  ^  Re{kx] »  we  can  quantify 
the  relationship  between  buoyancy  b  and  mh  in  (Bl) 


—iN^auw _  -A^g 

i(V<.  -  Ke)/:x2  '  f 

/+  2 


r(Vg  - 

L  If 


+  l\U\\ 


(B8) 


independent  of  the  direction  of  across-slope  propagation,  so  the  phase 
relation  between  U\\  and  b  becomes 


=  Arctan(  -) 

(y  e~  ^e)kz^ 


(B9) 


for  topographic  waves  that  are  either  propagating  and  standing  across  slope. 
The  phase  relation  between  mh  and  b  is  perturbed  from  the  inviscid  solution 
by  half  the  amount  that  viscous  near-inertial  waves  are  perturbed  from  their 
inviscid  solution.  We  apply  viscous  phase  relation  offsets  (A6)  and  (B9)  to 
seamount-  and  vortex-trapped  waves  under  the  assumption  that  seamount 
axisymetric  geometry  is  not  a  major  consideration. 
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